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Abstract 

We approach the analysis of dynamical and geometrical properties of nonholonomic mechan- 
ical systems from the discussion of a more general class of auxiliary constrained Hamiltonian 
systems. The latter is constructed in a manner that it comprises the mechanical system as 
a dynamical subsystem, which is confined to an invariant manifold. In certain aspects, the 
embedding system can be more easily analyzed than the mechanical system. We discuss 
the geometry and topology of the critical set of either system in the generic case, and prove 
results closely related to the strong Morse-Bott, and Conley-Zehnder inequalities. Further- 
more, we consider qualitative issues about the stability of motion in the vicinity of the critical 
set. Relations to sub-Riemannian geometry are pointed out, and possible implications of our 
results for engineering problems are sketched. 



This is a shortened and modified version of the author's PhD thesis in mechanical engineering at ETH 
Ziirich, from January 1999 (ETH-Diss 13017). The work was carried out at the institute of mechanics (now 
IMES), ETH Ziirich, before October 1998. 



Introduction 



This text is a shortened and modified version of the author's PhD thesis in mechanical 
engineering at ETH Ziirich, (ETH-Diss 13017, January 1999). Its subject matter is located in 
the area of theoretical mechanics, and discusses a particular class of constrained Hamiltonian 
systems, which contains the subclass of nonholonomic mechanical systems. We will be mainly 
concerned with geometrical and global topological issues about the critical manifolds, the 
stability of equilibria, aspects about numerical stability, and technical applications. 

An area of the engineering sciences, where nonholonomic mechanical systems are typ- 
ically encountered, is the field of multibody systems simulation. There is a multitude of 
different approaches to the description and analysis of such systems, stemming from various 
subareas of application. The geometrical approach given here is influenced by the work of 
R. Weber on Hamiltonian systems with nonholonomic constraints, and very essentially 
by the formalism developed by H. Brauchli based on an orthoprojector calculus, [0, |T2|, [3T[] . 
A construction for the Lagrangian case, which is closely related to what will be presented in 
section 3, has been given by Cardin and Favretti, |T^. A different approach in the Hamil- 
tonian picture is dealt with in the work of J. Van der Schaft and B.M. Maschke, A 



geometrical theory of nonholonomic systems with a strong influence of network theory has 
been developed by H. Yoshimura and T. Kawase, [0. The geometrical structure of non- 
holonomic systems with symmetries and the associated reduction theory, as well as aspects 
of their stability theory has been at the focus in the important works of J.M. Marsden and 
his collaborators, see for instance ^ 

In section one, we introduce a particular class of constrained Hamiltonian systems. 
Given a symplectic manifold (M, u) and a symplectic distribution V, we focus on the flow 
$t generated by the component of the Hamiltonian vector fleld Xh in V. 

Section two addresses the geometry and global topology of the critical set C of the 
constrained Hamiltonian system. The main technical tool used for this purpose is a gradient- 
like flow (pt, whose critical set € is identical to that of $j. Assuming that the Hamiltonian 
if : M — > R is a Morse function, it is proved that generically, C is a normal hyperbolic 
submanifold of M. The generalization of Morse theory developed by C. Conley and E. 
Zehnder is used to prove a topological formula for closed, compact that is closely related 
to the Morse-Bott inequalities. A second proof is given, based on the use of the Morse- Witten 
complex in Morse theory, to further clarify certain structural aspects. 

In section three, we address questions about the stability of the constrained Hamilto- 
nian systems in discussion, and conjecture a stability criterion for the marginally stable case. 
Relations to sub-Riemannian geometry are pointed out. 

In section four, we introduce Hamiltonian mechanical systems with Pfaffian constraints, 
which represent the case of highest technical relevance. It is shown that the physical orbits 
of such a system are conflned to a smooth submanifold Tlphys C T*Q. We construct an aux- 
iliary constrained Hamiltonian system of the type introduced in section one, which exhibits 
^phys as an invariant manifold, on which it generates the physical flow. Furthermore, we 
study the global topology of the critical manifold of the physical constrained system, and 
analyze the stability of equilibria of such systems. Finally, we propose a method to numeri- 
cally construct the generic connectivity components of the critical manifold. 
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I A NONINTEGRABLE GENERALIZATION OF 
DIRAC CONSTRAINTS 

Hamiltonian dynamics. Let {M,uj,H) be a Hamiltonian system. M denotes a smooth, 
symplectic 2n- manifold witli a symplectic structure uj G A^(M). H G C°^(M) is the 
Hamilton function, or the Hamiltonian in brief. For p = 1, . . . , 2n, A^(M) is the C°°{M)- 
module of p- forms on M. The Hamiltonian vector field Xh G r(TM) is determined by 

ixjj = —dH , 

where i stands for interior multiplication by a vector field. Given a smooth distribution 
W C TM, r{W) will denote the C°°(M)-module of smooth sections of W. The Hamiltonian 
flow is the 1-parameter group $t G Diff(M) generated by X^, with t G R, and $o = id- Its 
orbits obey 

dMx) = XniMx)) (1) 

for all X G M, and t G R. 

iu induces a smooth, non-degenerate Poisson structure on M, the R-bilinear, antisym- 
metric pairing on C°°(M) given by 

{f,g} = cu{Xj,Xg) . (2) 

It is a derivative in both of its arguments, and satisfies the Jacobi identity 

{/, {9,h}} + {h,{f,g}} + {g,{hJ}} = 0, 

thus (C°^(M), {■,■}) is a Lie algebra. Along every orbit of the Hamiltonian flow, 

dJiMx)) = {HJ}{^,{x)) (3) 

is satisfied for all / G C°°(M), and all x G M, t G R. 

Dirac constraints. Let j : M' — ^ M be an embedded, smooth, 2 /c- dimensional symplectic 
submanifold of M, endowed with the pullback symplectic structure j*uj. The induced Poisson 
bracket on M' is defined by 

{f,g}n = ifiu)iXf,Xg) 

for any pair of extensions f,g& C°°{M) of /, (7 G C°°{M'), and is in this context referred to 
as the Dirac bracket. 

If M' is locally characterized as the locus of common zeros of some family of functions 
Gi G C°^(M), with i = 1, . . . , 2(n — fc), it is possible to give the following exphcit construction 
of the Dirac bracket, [El]. Because M' C M is symplectic, the (n — k)"^ quantities 

Dij := {Gi,Gj} 

define a smooth, matrix- valued function that is invertible everywhere on M'. The explicit 
formula for the Dirac bracket is given by 

{f,9}D = {1,9} - {f,G,}D^^ {G,r9}, (4) 
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where D*-' denotes the components of the inverse of [D 



A natural generalization. We will next mimick this construction in the following more 
general setup. 

A distribution V over the base manifold M will be called symplectic if \4 is a sym- 
plectic subspace of T^M for all x G M with respect to uj^- Accordingly, V-^ denotes the 
distribution characterized by the property that its fibres are the symplectic complements of 
those of V with respect to cu, and is symplectic. Smoothness of V and u implies smoothness 
of V-^. If V is symplectic, the Whitney sum bundle V © V-^ equals TM. 

Thus, let V denote an integrable, smooth, symplectic rank 2A;-distribution V over M. 
Then, any section X G T{TM) can be written as 

X = X^ + X^" , 

where X are sections in r(\/(^)). Moreover, u{X^,X^^) = 0, by definition of V^. 
Associated to this decomposition of vector fields into components lying in V and V-'-, there 
exists an cj-orthogonal tensor TTy '■ TM — > TM with 

Ker(7ry) = , 7rv{X) =X VA G T{V) , 

which satisfies 

uj{7rv{X),Y)=uj{X,7Tv{Y)) 

for all A, y G r{TM). It will be referred to as the tu-orthogonal projector associated to V. 

Given a local family of spanning vector fields Yi, . . . ,Y2k for V, one can obtain an 
explicit formula for vry in a construction similar to ( To this end, one uses the fact 
that symplecticness of V is equivalent to the condition of invertibility of the matrix Cij := 

cuiY,,Y,). 

Lemma I.l Let Cij := uj{Yi,Yj), and let C^^ denote the components of its inverse. Define 
the 1-forms 6j{-) '■= i^(Yj, ■). Then, locally, ttv = C^^Yi ® 6j. 

Proof. Clearly, the range of C^^Yi 6j is V. Furthermore, 

{C'Yi ® = C'Y, ® 

because of O-iiYj) = Cij, and CijC^^ = 5'. To prove cu-orthogonality for the asserted expres- 
sion for TTy, consider an arbitrary pair of vector fields A, A', and observe that 

a;(7r^(A),A') = uj{C%iX)Yk, X') 
= C%iX)LuiYk,X') 

= c''eiix)ekix') 
= -c^''ei{x)e,{x') 

= -uj{Yi,X)&%{X') 

= uj{x,c"'ek{x')Yi) 

= UjiX,7TviX')), 

which establishes the claim. ■ 
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As a generalization of the usual Hamiltonian flow, the following dynamical system can 
be naturally associated to (M, a;, if, V). The component of Xh in V, 

Xl ■= MXh) C ViV) , 

generates a 1-parameter group of diffeomorphisms $4. The orbits of obey 

dA{x) = X^i^tix)) , (5) 

and are tangent to V for all initial conditions x e M. Embeddings R — M that are tangent 
to V are called V^-horizontal. If the condition of integrability imposed on V is dropped, this 
dynamical system will allow for the description of nonholonomic mechanics. 

Integrability of V implies that M is foliated in terms of 2/c-dimensional integral man- 
ifolds of V, which are, due to the symplecticness of V, symplectic submanifolds of M. 
Therefore, every leaf j : M' — M is an embedded invariant manifold under the action of $t. 

The induced dynamical system on any fixed leaf is equivalent to the puUback Hamil- 
tonian system {M, j*u!, H o j) that has been considered before in the discussion of Dirac 
constraints. In this sense, ( ^) generalizes the notion of Dirac constraints. 

Non-integrability. A new class of dynamical systems is obtained by discarding the re- 
quirement of integrability on V. 

By definition, the distribution V is non-integrable if there exists a filtration 

C C 1^2 C ■ ■ ■ C K, (6) 

inductively defined hj Vq = V, and 

V, = V,.^ + [Vo,V,.,], 

which is non-trivial, that is, Vi 7^ Vq. [■, ■] denotes the Lie bracket, and the sequence {Vi}{ 
is called the flag of V. For a pair of smooth distributions W, W, the distribution [W, W] is 
locally given by the linear span of all Lie brackets that can be taken between sections of W 
and W. If the fibre ranks of all Vi are base point independent, V is called equiregular. The 
smallest number r{V) at which the flag of V stabilizes, that is, for which 

K = K(y) \fs>r{V), 

is called the degree of non-holonomy of V. If Vr(v) = TM, V is said to satisfy Chow's 
condition, or to be totally non-holonomic. 

Proposition I.l (Frobenius condition) The symplectic distribution V is integrable if and 
only if in local coordinates, 

4- ■■= {drircy f, - d^ircy)';.) = (7) 

holds everywhere in M. 

Proof. Clearly, V is integrable if and only if vfy ([7ry(X), 7rv'(y)]) = is satisfied by all 
sections X, Y of TM, which is equivalent to Vi = V. The asserted formula is the local 
coordinate representation of this condition. ■ 
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Local coordinate representation. In a local Darboux chart with coordinates x*, the 
symplectic structure is represented by the matrix 



J 




-1. 



(8) 



Writing x{f) for the coordinate components of $t(a:o) with xq G M, the equations of motion 
( H) are represented by 



dtx{t) = (P, J d,H) {x{t)) 
= [J Pi d^^H) {xit)) 



(9) 



P denotes the matrix of vry, and is its transpose. Due to lemma |]T], P can be explicitly 
determined once one picks a local spanning family of vector fields for V . 



Almost Kahler structure. Let g denote a Riemannian metric on M. An almost complex 
structure J is a smooth bundle isomorphism J : TM TM with = —1. Together with 
g, it defines a two form satisfying 

Ug,j{X,Y) = giJX,Y) (10) 

for all sections X, Y E r(TM). The Riemannian metric g is hermitean if g{JX, JY) = 
g{X,Y), and Kdhler if Ug^j is closed. Thus, if Ugj = u (the symplectic structure), g is 
Kahler. A triple [g, J,uj) of this type is called compatible. 

Proposition 1.2 There is a compatible triple {g, J,uj), where u is the symplectic structure, 
such that P is symmetric with respect to g, and iTyJX = JiryX for all X G r(TM). 

Proof. Let us pick a smooth Riemannian metric g on M, relative to which P shall be 
symmetric. For instance, one may choose an arbitrary Riemannian metric g' on M, and 
use it to define g{X,Y) = g'{nvX,7ivY) + g'{7ivX,7tvY), where vfy denotes the projection 
complementary to vry. 

Since u is non-degenerate and smooth, there is a non-degenerate, smooth bundle map- 
ping K which is defined by 

u;{X,Y) = ~g{KX,Y). (11) 

It is skew symmetric with respect to g, because the symplectic structure is antisymmetric, 
so the adjoint operator of K with respect to g is K* = —K. 

The bundle map K can be used to construct an almost symplectic structure. Since 
K*K = —K^ is smooth, positive definite, non-degenerate and symmetric, there is a unique 
smooth, positive definite, symmetric bundle map A defined by = —K'^, which commutes 
with K. This immediately implies that the bundle map J = KA^^ satisfies = —1, so 
that it defines an almost complex structure. This is the standard proof for the fact that 
every almost symplectic manifold admits an almost complex structure, see for instance ||13|| . 

Because A is positive definite and symmetric, one can define a new metric by g{X, Y) = 
g{AX,Y), which obviously satisfies 

u;{X,Y)=g{JX,Y). (12) 
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Moreover, this metric is hermitean, since 



g{JX,JY) = ~g{KX,A-'KY) 
= -~g{X,K^A-'Y) 

= 9{X,Y). 

Because u is closed, g is not only hermitean, but even Kahler. 

The projector vry is symmetric with respect to g, as one concludes from the following 
consideration. The fact that P commutes with K follows from 

g{KnvX, Y) = cj(7ryX, Y) = uj{X, TCyY) = g{KX, TTyF) = g{7fvKX, Y) 

for all X and Y, since vry is symmetric with respect to g. P commutes with K, so also 
commutes with K^; thus it commutes with A"^. The linear operator A"^ is positive definite 
and symmetric, therefore commutativity with the symmetric operator P implies that P also 
commutes with A. It immediately follows that P is symmetric with respect to g. 

Finally, let us prove that vry commutes with J on r(TM). To this end, we observe 

that 

giJnvX, Y) = uj{ttvX, Y) = uj{X, vryF) = g{JX, TTyY) = ^(vry JX, Y) 

holds for all X and Y, as a result of which one concludes that nyJX = JiTyX is satisfied 
for all X, as claimed. This implies that J restricts to a bundle map J : V V. ■ 

If the J is covariantly constant with respect to the Levi-Civita connection associated 
to g, it is a complex structure. If uj admits a compatible Kahler metric g together with a 
complex structure J, (M, g, J) is a Kahler manifold. 

Symmetries. Given a Lie group G, {M,uj,H) admits a symplectic G-action \E' : G — > 
Diff(M), if 

\I^^u; = u; , H = H 

is satisfied for all h E G. Accordingly, for the case of a constrained Hamiltonian system, we 
will say that (M, u,H,V) possesses a G-symmetry if '^h*V = V holds for all h E G. 

Generalized Dirac bracket. In the same spirit in which one defines the Poisson structure 
induced by u, it is possible to introduce a smooth, R-bilinear, antisymmetric pairing on 
G°°(M) that is associated to the triple {M,u, V). In fact 

{f,g}v := u;{nv{Xf),nv{Xg)) (13) 

is a straightforward generalization of the Dirac and Poisson brackets. Along orbits of <I>t, 
one has 

for all X G M, in analogy to ( ^). The bracket ( [13|) does not satisfy the Jacobi identity if V 
is non-integrable, but it satisfies a Jacobi identity on every (symplectic) integral manifold if 
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V is integrable. 

Energy conservation and symplecticness. Let us finally mention two further important 
properties of tfie flow $t G DifT(M), t G R, of the constrained system. 

Proposition 1.3 The energy H is an integral of motion of the dynamical system (\^. 

Proof. This follows from the antisymmetry of the generalized Dirac bracket, which implies 
that dtH = {H, H}v = 0. ■ 

Finally, let us consider 

= [dixvu + ix^duj^ 

= -^;d{ {Pv)l diH dx^ ) 

= - {di{Pv)idiH dx^ Adx'') 

= -^t (^l {di{Pv)l-dk{Pv)i) diHdx'Adx'^ . 
This result implies that the restriction of dt^^u to X, F G r(V") satisfies 

dt^iuix, Y) = (I di, m x'Y' 
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where is defined in lemma |]T|, which states that the the right hand side vanishes 
identically if and only if V is integrable. 

Thus, if V is integrable, the restriction of to x V" equals its value for t = 0, given 
by ct;(7ry(-), 7ry(-)). Hence, on every integral manifold j : M' ^ M oiV, is symplectic 
with respect to the puUback symplectic structure j*uo. 
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II THE GEOMETRY AND TOPOLOGY OF THE 
CRITICAL MANIFOLD 

The critical set C of a dynamical system is defined as the set of equilibrium solutions. The 
main purpose of this section is to study geometrical and global topological properties of the 
critical sets €. exhibited by constrained Hamiltonian systems of the type {M,u), H,V). An 
application of Sard's theorem will demonstrate that generically, £ is a smooth 2{n — k)- 
dimensional submanifold €.gen C M, the critical manifold. 

Moreover, it will be demonstrated that there exists a precise correspondence between 
€gen and the critical points of if if : M — > R is picked to be a Morse function. In 
particular, we will prove that for compact M without boundary, the Poincare polynomials 
of M and those of the connectivity components of ^gen are related in a way reminiscent of 
the Morse-Bott inequalities. 

Two different proofs will in fact be presented. The first proof is based on the gener- 
alization of Morse and Morse-Smale theory ^ developed by C. Conley and E. 



Zehnder in |T^, applied to an auxiliary gradient-like system. The second proof is based on 
the comparison of the Morse- Witten complexes of {M,H) and {(tgen, Hltr^^^). 

The special case of mechanical systems (where M is noncompact) will be analyzed in 
the last section. 

The following standard definitions are necessary for the subsequent discussion. 

Definition II. 1 The zeros of the one form dH are called critical points of H, and the value 
of H at a critical point is called a critical value. A level surface T^e that contains a critical 
value E of H is called a critical level surface. A level surface S e that contains no critical 
points of H is called regular, and consequently, E is then called a regular value. 

II. 1 Generic properties of the critical set 

The critical set of the constrained Hamiltonian system (M, u;, if, V) is given by 

C = {xeM\ X^{x) = 0} C M . 

The following theorem holds independently of the fact whether V is integrable or not. 

Theorem II. 1 Generically, the critical set is a piecewise smooth, 2{n — k)- dimensional 
submanifold of M. 

Proof. Let denote a smooth, local family of spanning vector fields for V over an 

open neighborhood U C M. Since V is symplectic, the fact that Xj^ is a section of V 
implies that uj{Yi,X'^) cannot be identically zero for all i and everywhere in U. Due to the 
cu-orthogonality of vry, and vryFj = Yi, 

uj{Y,,X];) = uj{ttv{Y^),Xh) = uj{Y,,Xh) = Y,{H) . 

Thus, with F := {Yi{H), Y^kiH)) G C°°(f/, R^'^), it is clear that Cnt/ = F-\Q). Since F 
is smooth, Sard's theorem implies that regular values, having smooth, 2(n — /i;)-dimensional 
submanifolds of U as preimages, are dense in FjU) ■ 
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Next, let us pick a local spanning family {Yi E T(y)}f^-^ for V that satisfies 



with J :- 



. This choice is always possible. 



Ik 
-Ik 

Furthermore, introducing the associated family of 1-forms {9i} by 9i{-) := Lj{Yi,-), 
theorem ^ implies that 

Try = T% ® Oj , 

where J^^ are the components of = — J. Expanding with respect to the basis {Yi} 
gives 



X 



H 



TTviX^) = - Y,{H) Y, , 
where one uses the relationship 6j{Xj^) = Yj{H) obtained in the proof of theorem |IL1 



(14) 



Proposition II. 1 Under the genericity assumption of theorem the 2k x 2k-matrix 

[Yis{Yi{H)){xo)] is invertible for all Xq G and every local spanning family {Yi G r(V)} of 
V. 

Proof. Let us pick a local basis for V, and {Zj}^^^~^^ for V-^, which together span 



TM. Let xo G Cgen, and assume the generic situation of theorem p.l| . Because Cgen is 
defined as the set of zeros of the vector field ( |T^, the kernel of the linear map 



dF^nr'^Y 



where Fi := Yi{H), is precisely Tx^^gen, and has a dimension 2{n — k). 

In the basis {FiUo, • • • > ^2fcU'o, ^lUo, • • • , ^2(n-fc)Uo}, its matrix is given by 



A 



[Av A, 



where Ay := [Yi{Fj) J^^Ykl^^], a.nd Ay ± := [Zi{Fj)J^^Yk\^g]. Bringing A into upper triangular 
form. Ay is also transformed into upper triangular form. Because the rank of A is 2k, 
and Ay is a 2A; x 2A;-matrix, its upper triangular form has 2k nonzero diagonal elements. 
Consequently, Ay is invertible, and due to the invertibility of J, one arrives at the assertion. 



Corollary II. 1 Let €.gen satisfy the genericity assumption of theorem |//. 1\ . If V is inte- 
grable, the intersection of any integral manifold of V with (tgen is a discrete set. 

Proof. The previous proposition implies that generically, integral manifolds of V intersect 
<tgen transversely. Their dimensions are complementary, hence the intersection set is zero- 
dimensional. ■ 
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II. 2 Topology of the critical manifold 



The Hessian of the constrained system. The usual coordinate invariant definition of 



the Hessian of H corresponds to VdH, evaluated at its critical points |]22[, where V denotes 
the Levi-Civita connection of the Kahler metric g. 

A generalized Hessian for the constrained Hamiltonian system can easily be defined 
along these lines. Let us write vTy for the dual projector associated to vry, which acts 
on sections of the cotangent bundle T*M. Let 6 denote any one-form. Then of course, 
{7iy9,X) = {OjTTvX). In any local chart, the matrix of VTy is the transpose of the matrix of 

TTy. 

The obvious generalization of the Hessian is the tensor V {TCydH) . It acts as a bilinear 
form on pairs of vector fields in terms of 

V{7rUH){X,Y) = {Vx{7rUH),Y) 

= ii{ny)iH,)^,-TUiny)iH,)X^Y% 

where F^^ are the Christoffel symbols. The second term in the bracket on the lower line 
is zero on C, because {7iv)iHj = on C The non-vanishing term in V^HydH) on £ is 
determined by the matrix 

Krs = {{TrvYrH,,),,. (15) 

A straightforward calculation shows that {nvYiKjk = Kik holds everywhere on C Obviously, 
the rank of K is bounded from above by the rank of Try, corresponding to the rank 2k of V . 

The corank of K\a equals the dimension of the connectivity component of C which 
contains a. This is because for any a' G C close to a, Taylor's theorem implies that 

{{'Kv)iH^,){a') = {{Tiy)lH,,){a) + K\a{a' - a) + 0{\\a! - af). 

The left hand side and the first term on the right hand are both zero, since €. is so defined, 
thus — s> in the limit ||a' — a|| — > 0. That is, the tangent space T^C is equal to keri^l^, 
hence the corank of K\a is indeed the dimension of the component Cj. 



II. 2.1 An auxiliary gradient-like dynamical system 

The physical fiow of the constrained Hamiltonian system turns out to of very limited use for 
the study of the global topology of This is because invariant sets of the physical fiow not 
only contain fixed points, but also periodic orbits. 

However, the auxiliary dynamical system on M defined by 

%(t) = -(7ryV,//)(7(t)), (16) 
for 7 : / C R — » M, is an extremely powerful tool for this purpose. Let us denote its fiow by 

0t G Diff(M) . 

Its orbits are of course tangent to V , and it is clear that $j and (pt share the same critical 
set £. 
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Definition II. 2 A flow is gradient-like if there exists a function / : M R that decreases 
strictly along all of its non-constant orbits. 



Proposition II. 2 The flow (pt is gradient-like. 

Proof. One can easily check that H decreases strictly along the non-constant orbits of 

m{^{t)) = { dH{^{t)) , dMt) ) 

= -giV.H, vrvV.if) (7(t)) 
= -ginvVgH, 7Tv\/gH){^{t)) 
< 0. 

We have here used the fact that {g, J,uj,ttv) is a compatible quadruple. This is the main 
motivation for the introduction of these quantities. ■ 

It is immediately clear that (pt generates no periodic trajectories, hence €. comprises all 
invariant sets of (pf 

The gradient-like fiow close to €. Let us briefly consider the orbits of (pt in a tubular 
e-neighborhood of (with respect to the Riemannian distance induced by the Kahler metric 
g). We pick an arbitrary element a G C and local coordinates x*, with the origin at a. The 
equations of motion ( 0) are given by 

dtx'-it) = - K:{a)x\t) + 0(||xf ), (17) 

where Kjs is the 'constrained Hessian' at a, and = g^^Kjs, and || ■ || is with respect to g 
at a. 

Pa will from now on denote the matrix of 7Ty{a). Furthermore, let Ip^ : Va ^ TaM 
denote the embedding. Moreover, we will write Aa = [Kl{a)] for the Jacobian matrix of 
TTyVg-ff at a in the present chart. From the ^f-orthogonality of ttv follows straightforwardly 

that PaAa = Aa- 

The linearized equations of motion at a read dtx{t) = —Aax{t), and satisfy the lin- 
earized constraints dtx{t) = Padtx{t). 

Let Q denote the ^f-orthoprojector onto kerA^ = TaC, and let Iq^ : Ta€. ^ TaM be the 
embedding. Its complement Qa = ^ — Qa projects T^M orthogonally onto the fibre Na^ of 
the normal bundle of £, and likewise, let Iq^ : TaN ^ TaM be the embedding. The orbits 
of the linearized system are then given by 

x{t) = exp{-tIp^AalQ^Qa)xo (18) 

= exp{-tIp^AalQ,Qa)lQaQaXo + Ig^QaXo, (19) 

where the initial condition is determined by a;(0) = xq. The vector QaXo G Ta€. approxi- 
mately connects a with some critical point close to a, while QaXo lies in Na€.. 
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II. 2. 2 Morse functions and non-degenerate critical manifolds 

Let us recall some standard definitions from Morse- and Morse-Bott theory that will be 
needed in the subsequent discussion. 

Definition II. 3 The dimensions of the zero and negative eigenspaces of the Hessian of f 
at a critical point a are called the nullity and the index of the critical point a. If all critical 
points of f : M H have a zero nullity, f is called a Morse function, and the index is then 
called the Morse index of a. 

If the critical points of f are not isolated, but elements of critical manifolds that are 
non-degenerate in the sense of Bott, H is called a Morse-Bott function 

Throughout this section, we will assume that is a Morse function. 

Definition II. 4 A connectivity component is locally normal hyperbolic at the point a G C 
with respect to (pt if it is a manifold at a, and if the restriction of Aa to the normal space Na€. 
is non- degenerate. A connectivity component di is called non-degenerate if it is a manifold 
that is everywhere normal hyperbolic with respect to (f)f The index of a non- degenerate 
connectivity component is the number of eigenvalues of the constrained Hessian Aa on €,i 
that are contained in the negative half plane. 



Proposition II. 3 In the generic case, the critical manifold is normal hyperbolic with respect 
to the gradient-like flow. 



Proof. This follows straightforwardly from lemma ILl 



Let Cj, i = 1, ..,1 denote the connectivity components of C = UCj, and let ji : (ti ^ M 
denote the embedding. 

Proposition II. 4 Assume that €. satisfies the genericity assumption in the sense of Sard's 
theorem. Then, if H : M ^ H is a Morse function, := H o : — > R zs also a Morse 
function. A point x E €i is a critical point of if and only if it is a critical point of H. 

Proof. The fact that every critical point of is a zero locus of TTyVgi^la, and thus an 
element of <t, is trivial. 

To prove the converse statement, let £j be generic, so it has a dimension 2{n — k). 
Normal hyperbolicity implies that the restriction of Palq^QalpaPa to the normal space Na^i 
defines an invertible map for all a G Cj. Let a now denote an extremum of H\€,i. Then, 

{dH,v)l = giygH,v)l = Q 

holds for all v G TqCj. Consequently, V gH\a is a vector in Na^, and thus, V gH\a = QaVgH\a. 
Moreover, by definition of £, the condition PaVgH\a = is satisfied, which in turn implies 
QalpaPalQaQaV gH\a = 0. Heucc, siucc Qalp^PalQaQa ■ Na^i ^ Na^i is iuvertiblc, it follows 
that S/gH\a = 0. 

The Hessian of the restriction of H at any critical point of Hi is always nondegenerate, 
thus Hi is by itself a Morse function on CCj. ■ 
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Corollary II. 2 The critical points of H\(rg^^ are precisely the critical points of H : M — > R. 
Let €j he a non- generic connectivity components with manifold structure that is normal 
hyperbolic. Then, C T,H{€i)- 

Proof. The first assertion follows trivially from the previous proposition. 

Next, let us assume that £j is some non-generic connectivity component of £ with 
manifold structure that is normal hyperbolic (which is in no sense necessitated by the sys- 
tem). Then, the proof of the previous proposition can be applied to the situation for CCj, 
and demonstrates that there are no extrema of H\it^. Thus, is a submanifold of the level 
surface S//(ev). ■ 



II. 2.3 Index pairs and relative (co-) homology 

Let us next focus on the global topology of for assumptions that are weaker than genericity. 
In fact, we will only assume that the generic connectivity components Cgen C are non- 
degenerate, compact and without boundary, and that €\ C^en is compact. In this situation, 
a generalization of the Morse-Bott inequalities can be derived for based on an application 
of the theory of C. Conley and E. Zehnder for flows on Banach spaces with compact invariant 
sets. We will only give a very short survey of Conley-Zehnder (CZ) theory, and refer the 



reader to the original article [ITS 



Definition II. 5 Let he any compact component of <t. An index pair associated to €i is 
a pair of compact sets {Ni,Ni) that possesses the following properties. The interior of Ni 
contains €i, and moreover, CCj is the maximal invariant set under 0t iiT' the interior of Ni. 
Ni is a compact suhset of Ni that has empty intersection with <ti, and the trajectories of all 
points in Ni that leave Ni at some time under the gradient-like flow (pt intersect Ni. Ni is 
called the exit set of Ni. 

The p-th relative homology group Hp{X,A) of a pair of manifolds A G X is defined 
as follows. Let Cp{X) be the group of p-cycles of X. Since Cp{A) is a subgroup of Cp{X), 
the quotient group Cp{X)/Cp{A) is well-defined. Let Bp{X,A) denote the subgroup of 
Cp{X) I Cp(y4) which consists of boundaries. The associated quotient group is the j9-th relative 
homology group Hp{X, A). 

The p-th relative de Rham cohomology group Hp{X, A) consists of cohomology classes 
that are represented by closed p-forms on X whose restriction to A (via the puUback of 
the embedding A — > X) is exact. If X and A are both orientable, the coefficients of the 
relative homology and cohomology groups can be picked from R or Z, and otherwise from Z2. 

The relative cohomology of an index pair. It was proved in [|1^ that the homotopy 
type of the pointed space Ni/Ni only depends on £j, so that H*{Ni,Ni) is independent 
of the particular choice of index pairs (the space Ni/Ni is obtained from collapsing the 
subspace Ni of Ni to a point). The equivalence class [Ni/Ni] of pointed topological spaces 
under homotopy only depends on and is called the Conley index of In the present 
analysis, we will only consider flows exhibiting normal hyperbolic critical manifolds. In this 
special case, the result essentially reduces to a generalization of Morse-Bott theory. We will 
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now determine the relative cohomology H*{Ni, Ni) of an index pair [Ni, Ni) associated to a 
generic connectivity component <ti. 

Proposition II. 5 Let C €gen be a generic connectivity component that is compact and 
without boundary, and let {Ni,Ni) denote any associated index pair. Then, 

H'^+^^{N,,N,) = H^idi) , (20) 

where fii is the index of Ci, and q = 0, . . . , dim{(ti) . 

Proof. Let us consider, for some sufficiently small eo > 0, a compact tubular eo-neighborhood 
U of <ti (of dimension 2n), and let 

W[r{iti) := {W'i€i)U€i)nU 

denote the intersection of the center unstable manifold of <ti with U. W~{€i) denotes the 
unstable manifold of Cj. Pick some small, positive e < eo, and let 11^ be the compact tubular 
e-neighborhood of Wij^{Ci) in U. 

It is clear that letting e continuously go to zero, a homotopy equivalence of tubular 
neighborhoods is obtained, for which W[j^{€i) is a deformation retract. Let 

ur' ■■= du, n 0R(f/,) 

denote the intersection of dU^ with all orbits of the gradient-like flow that contain points in 
U<^. Then, evidently, {U^,U°^^) is an index pair for £j, and by letting e continuously go to 
zero, is homotopically retracted to dW[j^{<ti). 

Thus, homotopy invariance implies that the relative cohomology groups obey 

H*{u,,ur') = /j*(iy^"(£,),9iy^"(c:,)). 

Due to the normal hyperbolicity of CCj with respect to the gradient-like flow, W^^{<^i) has a 
constant dimension Ui + /i(£i) everywhere, where nj = dimCj. Therefore, one obtains from 
Lefschetz duality [|I^ that 

where /ij = fi{^i), the index of £j. It is clear that £j is a deformation retract of the interior 
of Wff^{Ci), so that the respective cohomology groups are isomorphic. 
Due to dim{(ti) = n^, Poincare duality implies 

so that with q := Ui — p, 

if^+^^(f/e, t/r*) = H''{(t,), (21) 

which proves the claim. ■ 
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II. 2.4 The Conley-Zehnder inequalities 

Here follows a very rough summary of the main steps in the derivation of the Conley-Zehnder 
inequalities, [|l^ . 



Definition II. 6 Let I denote a compact invariant invariant set under (pf A Morse de- 
composition of I is a finite, disjoint family of compact, invariant subsets {Mi, . . . , M„} that 
satisfies the following requirement on the ordering. For every x G / \ UjMj, there exists a 
pair of indices i < j, such that 



lim (f)t{x) C Mi lim (^t(x) C Mj. 

t— >— oo t^ca 



Such an ordering, if it exists, is called admissible, and the Mi are called Morse sets of I . 



A central result proved in is that for every compact invariant set I admitting an 
admissibly ordered Morse decomposition, there exists an increasing sequence of compact sets 
A^j with No C Ni C . . . C Nm, such that {Ni, A^j_i) is an index pair for Mi, and {Nm, Nq) is 
an index pair for /. 

Consider compact manifolds A D B D C . The exact sequence of relative cohomo logics 

. . H''{A, B) H''{A, C) H\B, C) ^ H^+\A, B) ^ . . . 

implies, in a standard fashion known from Morse theory, that, with r^^p denoting the rank of 
H^i^Ni, iVi_i), the identity 



holds for the indicated Poincare polynomials Here, Bj is the j-th Betti number of the 
index pair (iV^, Nq) of I, and Q(A) is a polynomial in A with non-negative integer coefficients. 
These are the strong Conley-Zehnder inequalities; a corollary that straightforwardly follows, 
due to the positivity of the coefficients of Q(A), is that 



Er- > B 



holds. These are the weak Conley-Zehnder inequalities. 

If the symplectic manifold M is compact and closed, and if £ is assumed to be non- 
degenerate, the following results arise from application of the CZ inequalities. The invariant 
set I can be chosen to be equal to M. So we let Nm = M and A^'o = denote the top 
and bottom elements of the sequence defined above. Furthermore, we order the connected 
elements of C according to the descending values of the maximum of H attained on each Cj. 
Then it follows that €. furnishes a Morse decomposition for M. The homology groups of M 
are isomorphic to the relative homology groups of the index pair {Nm, Nq). So the numbers 
Bp are the Betti numbers of the compact symplectic manifold M. 

The number rj p is the rank of the p-th relative cohomology group of the non-degenerate 
critical manifold C,, the index i being determined by the Morse decomposition. From ( pTf ), 
we deduce that 

ri,p = dimH'^P-'''{<ti) 
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(we recall that /Xj is the index of (ti). In other words, r^^p is the (p — /ii)-th Betti number 
of Cj, written Bi^p^^^ in brief. Assuming that the number of connected components of €. is 
finite, the Conley-Zehnder inequalities thus result in 

B^,p X^^^' = J2 Bp\'' + {l + \) Q(A), (22) 

i,p p 

which implies that 

Bi-p-f,, > Bp (23) 

i 

is satisfied. So the global topology of M enforces lower bounds on the ranks of the homology 
groups of the connected components of €. Setting the variable A equal to —1, one obtains 
that 

i,p i 

where x denotes the Euler characteristic. 

Remark. In the case of mechanical systems, the phase space of the relevant constrained 
Hamiltonian system is non-compact, and the critical manifold is generally unbounded. There- 
fore, the arguments used here do not apply. However, since in that case, M and €. are vector 
bundles, we are nevertheless able to prove results that are fully analogous to ( P^ . 

Generic connectivity components. One can actually prove a stronger result than ( 
given the special structure of the system at hand. 

Proposition II. 6 Assume that C\ C^en is a disjoint union of -manifolds. Then, 

J2 B.,,p A^'+^' = ^ A^ + (1 + A) Q(A) . (24) 

Bi^p are the p-th Betti numbers of the connectivity components of €gen, Bp are the Betti 
numbers of M, and Q is a polynomial with non-negative integer coefficients. 

Proof. We will show that can be suitably deformed such that the components of €\<tgen 
can be smoothed away. To this end, we construct an auxiliary, continuous vector field X^, 
corresponding to a small deformation of the gradient-like vector field TTyVgH. We pick a 
small positive, real number e, and consider the compact neighborhoods 

U,{€i) = {xe M \ distg{x, <ti) < e} (25) 

of connectivity components £j C £ \ Cgen- Here, distg denotes the Riemannian distance 
function induced by the Kahler metric g. The vector field X^ shall be given by nyVgH in 
M \ Ue{Ci), and in the interior of every f/e(Cj) with C \ Cgen, 

X,\, = 7ivVgH\, + eh{x)VgH\, . (26) 
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Here, we have introduced a C^-function h : U^i^i) — > [0, 1] obeying h\it^ = 1 and h\dUt{€i) = 0, 
that is strictly monotonic along the the flow lines generated by TTyVgH. Such a function h 
exists because C \ (tgen is a disjoint union of C^-manifolds. 

We have proved above that for all C € \ Cgen, ^ gH is strictly non-zero in U^{'ti). 
Now let us consider 

g{X,,VgH) = {\\7^yVgH\\l){x)+th{x){\\VgH\\l){x). 

Here we have used the ^f-symmetry of P, and \\X\\^g = g{X,X). The flrst term on the right 
hand side is non-zero on the boundary of Ue{€.i), while the second term vanishes. Moreover, 
the second term is non-zero everywhere in the interior of Ue{€i). Therefore, vanishes 
nowhere in Ue{€i). 

This implies that the vector fleld X^ is a deformation of iryVgH that only exhibits 
€gen as a critical set. Notice that the generic components cannot be removed in this manner, 
since they contain critical points of H. 

The scalar product of X^ with VgH is not only non-zero, but also strictly positive in 
every Ue{*ti), which shows that X^ is a gradient-like flow (the left hand side of the above 
expression corresponds to the time derivative of H along orbits of in U^{€.i)). 

Clearly, the order of magnitude of ||Xe — nvVgEWg is at most 0(e) everywhere on M. 
Consequently, it is possible to pick X^ arbitrarily close to TiyVgH in the || • ||oo-norm on 
r(TM) that is induced by || ■ \\g. 

Carrying out the Conley-Zehnder construction with respect to the flow generated by 
Xg yields ( |2^). This result does not require the assumption of normal hyperbolicity on (t. ■ 

II. 2. 5 Proof of ( ^) using the Morse- Witten complex 

We will now give a different proof of ( based on the existence of the Morse- Witten 
complex. Our motivation is, on the one hand, to clarify the orbit structure of the gradient- 
like system, and, on the other hand, to further comment on the "position" of <tgen in M 
relative to the equilibria of the free system (corresponding to the critical points of H). 

To this end, let us flrst make some brief general remarks about Morse theory, cf. 

In the basic setting, there is a Morse function / : M ^ R on a compact manifold M, 
without boundary, say. The index of an isolated critical point precisely is its Morse index, 
hence the strong Conley-Zehnder inequalities ( yield 

5^A*'Xp = 5^APi?p + (l + A)Q(A), 

p p 

where Np is the number of critical points of / with a Morse index p. Bp is the p-th Betti 
number of M, and Q{t) is a polynomial with non-negative integer coefficients. These are the 
strong Morse inequalities, from which the weak Morse inequalities Np > Bp are immediately 
inferred. The standard proof uses the gradient flow generated by —Vgf with respect to an 
arbitrary auxiliary metric. 

If / is not assumed to be a Morse, but more generally, a Morse-Bott function the 
gradient flow of the vector fleld — V^/ can again be used, but then to derive the strong 
Conley-Zehnder inequalities for the critical manifolds of /. 
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Morse theory, which has been a classical topic in differential topology since the 1930 's, 
experienced a tremendous new increase of activity in the early 1980's, due to the seminal work 
of E. Witten [^]. His proof of the weak Morse inequalities is based on a deformation and 



localization argument in Hodge theory, interpreted as supersymmetric quantum mechanics. 
Subsequently, the proof of the Morse-Bott inequalities using heat kernel methods has been 
given by J.-M. Bismut in 0. 

This new approach to the strong Morse inequalities is based on the differential complex 
constructed from the the set of critical points of /, which is often referred to as the Morse- 
Witten complex. At least Milnor, Smale and Thom independently have arrived at some 
form of it already earlier |T^, The (co-)homology of the Witten complex is isomorphic 
to the singular (de Rham co-)homology of M, and straightforwardly implies the strong 
Morse inequalities. A proof of the Morse and Morse-Bott inequalities that is based on the 
construction of the Morse- Witten complex is given in |^ . 

The Morse- Witten complex has been generalized to infinite dimensional systems by 



A. Floer WTl, which has led to extremely fruitful applications. For instance, the celebrated 



Arnol'd conjecture has been proved by use of Floer homology, A beautiful survey is 

given in |RB . 



We will now discuss an alternative proof of ( ^4]) that is based on the Morse- Witten 
complex for non-degenerate Morse functions. 

The Morse- Witten complex. Let M be a compact, closed, orientable and smooth man- 
ifold, together with a Morse function / : M — >■ R. We let denote the free Z-module 
generated by the critical points of / with a Morse index p. The set C = ©pC^ is the free 
Z-module generated by the critical points of /, which is graded by their Morse indices. There 
exists a natural coboundary operator 

5 -.CP ^ CP+^ 
that obeys 5^ = 0, which we will define below. 

Theorem II. 2 The cohomology of the differential complex (C,5) is isomorphic to the de 
Rham cohomology of M, ker5/im5 = H*{M, Z). 



A proof of this theorem has been given by Floer in [T7| based on Conley-Zehnder 



theory. Other proofs can be found in The argument in the original publication 

is based on the quantum mechanical tunneling effect. 



The coboundary operator is defined as follows ||, [T^, 0. We denote the unstable 
and stable manifold of a critical point a of f under the gradient flow by W~ and W^, 
respectively, and assign an arbitrary orientation to every W~ . Since M is assumed to be 
oriented, the orientation of W~ at every critical point a induces an orientation of W^. The set 
of Morse functions for which the stable and unstable manifolds intersect transversely is dense 
in C°°(M). Thus, we may generically assume that all W~ and intersect transversely. 
The dimension of W~ equals the Morse index /i(a) of a, and the dimension of the intersection 

M(a, a') ^ W~ n iy+ 

is given by max{/i(a) — /i(a'), 0}. 
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In order to define the coboundary operator, let us consider pairs of critical points a 
and a', wliose relative Morse index has the value 1, say fi{a) = p + 1 and /i(a') = p. It 
immediately follows that M{a, a') is a finite collection of gradient lines that connect a with 



a'. 



The intersection of M{a,a') with every regular level surface Sc of / with /(Sc) = c 
lying between /(a) and f{a') is transverse, and consists of a finite collection of isolated 



points. The hypersurface Sc, being a level surface of /, is orientable, so we pick the 
orientation, which, combined with the section V^/ of its normal bundle, shall agree with the 
orientation of M. 

In addition, the intersection both of W~ and with Sc is transverse, and the sub- 
manifolds 

^ w- n ^ n 

of Sc are smooth, compact and closed. In addition, their dimensions add up to the dimension 
of Sc. To every point b of the set M(a, a')nSc = W~^nW^ ^, one assigns the number 7(6) = 1 
if the induced orientation of 

agrees with the one picked for Sc, and 7(6) = — 1 otherwise. 
Definition II. 7 The sum 

{a,6a')= 7(6) 

beM{a,a')nT,c 

is the intersection number tl (W^cTc? ^J^c) of the oriented submanifolds and ^ of T^c, 
j^/ . If for a pair a and a' of critical points with a relative Morse index 1, this intersection 
number is nonzero, they will be said to be effectively connected (by gradient lines). 



Definition II. 8 The coboundary operator of the Morse-Witten complex is the 7i-linear map 
5:CP^ CP+i defined by 

6a' = (a, 5a') a. 

^{a)=p+l 

The coboundary operator satisfies 5"^ = 0, and its cohomology is isomorphic to the de 
Rham cohomology of M. Proofs of this statement can be found in |Q, |r^, The image 
of a critical point of Morse index p under the coboundary map consists of the critical points 
of Morse index p + 1 to which it is effectively connected. One also refers to the intersection 
number (a, 6a') as the matrix element of the 5-operator with respect to the elements a and 
a' of C. 

The existence of the Morse-Witten complex straightforwardly implies the strong Morse 
inequalities, as one deduces from the following fact. Let us denote the p-th cocyle group by 
2P c C^, which is defined as the intersection of ker5 with C^, and let B*^ C C^, the p-th 
coboundary group, denote the image of C^~^ under 6. Clearly, is a subset of because 
S is nilpotent, therefore the set = Z^\B^ is well-defined. It is the p-th cohomology group 



of the Morse-Witten complex. Since by theorem 11.2, the cohomology of the Morse-Witten 
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complex is isomorphic to the de Rham cohomology of M, the rank of coincides with the 
p-th Betti number Bp{M) of M. 

The image of in C^"^^ under 5 is given by B^+^. Denoting the preimage of B^+^ in 
by 5~^(B*'+^), which is isomorphic to B^+^, one has 

C^' = H^©BP©r^(Bf+^), 

so that 

dimC*" = Bp{M) + dimB*" + dimB*'+^ 

The dimension of equals the number Np of critical points of / with a Morse index p. 
Multiplying both sides of the equality sign with A^, and summing over p, one finds 

^A^ATp = J]A^'5p(M) + J]Af(dimB*' + dimBP+^) 

= ^AfSp(M) + (l + A)^A^'-MimBP (27) 

(notice here that both B° and B^"+^ are empty). These are the strong Morse inequalities, 
and the polynomial Q{t) defined at the beginning of this section now has a very definitive 
interpretation: 

Q(A) = ^A^'-MimBf. 

Evidently, dimB^' is the number of critical points of Morse index p that are effectively con- 
nected to critical points of Morse index p — 1 via gradient lines of /. 

Comparing the Morse- Witten complexes of {M,H) and (Cgen, ^^kgen)- Let us next 
relate the Morse- Witten complexes of {M,H) and {^gen, H\€gen) ^^ch other. To this 
end, it is necessary to recall that the generic connectivity components C^en of ^ contain all 
critical points of H ^ but no other conditional extrema, and that they are necessarily normal 
hyperbolic with respect to (f)f 

Let Ai := {oti,!, • • • , ai,m} denote the set of critical points of H that are contained in £j, 
and let ^{ai^r) be the associated Morse indices of if : M — > R. Furthermore, let Hi = H\(s^ 
denote the restriction of the Hamiltonian to (tj. We have previously shown that the map 
Hi : ^ R is a Morse function, whose critical points are precisely the elements of Ai. 
Furthermore, it is clear that the number of negative eigenvalues of the Hessian of H at any 
element of Ai, whose eigenspaces are normal to (Tj, equals n{Ci), the index of Cj. 

The Morse index of Ui^r with respect to Hi is thus ^{ai^r) — l^i^i)- The Morse- Witten 
complex associated to is defined in terms of the free Z-module generated by the elements 
of Ai, which is graded by the Morse indices p of the critical points of Hi, 

To define the coboundary operator 5i : C\ ^ C^"*"^, one uses the gradient flow on €j generated 
by Hi. One then concludes that 

keT5i/\-m5i^ H*{ti,Z). (28) 
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Application of ( shows that for every £j G €gen, 

J^A^iV,, = 5^A^'S,(C,) + (l + A)5^AP-MimBf, (29) 

p p p 

where is the p-th coboundary group of the Morse- Witten complex of (ti, and iVj^p is the 
number of critical points of Hi on €i whose Morse index is p. 

Since every critical point of H lies on precisely one generic component Cj, the number 
Nq of critical points of H with a Morse index q is given by 



^p = J2^-'P-f^ 



Multiplying both sides of ( ^9|) with A^^*^-*, and summing over i, one obtains 

J2 xt^i'^)+PN.^^= J2 A^('^)+P5p(Ci) + (1 + A) A'^('^)+P-MimB 



which becomes, after reindexing + P — ^ <?, 

^A''iV,= 5^ A^i?,_^(£,)(G:.) + (l + A) 5^ A^"MimB 

Combining this result with the strong Morse inequalities 

J2 X'Ng = J2 X'BgiM) + (1 + A) ^ A'^-MimB", 



one obtains that 



J2 \'Bq_,^^){<Li) = Y,>^^Bq{M) 
+ (1 + A) ^ X'^-^ I dimB"? - J2 diniB 



One observes that formula ( ^4|) implies that the polynomial which is multiplied by (1 + A) 
has non-negative integer coefficients. 

Conversely, if one can prove that for all q, 

dimB5> dimBr^^*^^ (30) 

holds, one would also obtain an alternative proof of ( |2^ ). The main observation here is that 
the left hand side is defined in terms of the operator 6 associated to (M, H), while the right 
hand side is defined in terms of the operators Si associated to all (£j. Hi) with G Cgen- 

The quantity dimB^"'^'^'^'* denotes the number of critical points of H with a Morse 
index q in CCj, which are effectively connected to critical points of Morse index p + 1 in CCj 
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via gradient lines of the Morse function Hi on Therefore, the sum on the right hand side 
of ( |3y) equals the number of those critical points of H with a Morse index q, which are 
effectively connected to critical points of Morse index g + 1 via gradient lines of the functions 
H o jj on all generic Cf, here, ji : €j M is the inclusion map. 

Proof of inequality ( |30|). We will now prove inequality ( ^) by relating the coboundary 
operators of the Morse- Witten complexes of {M,H) and (C^en, -f^kgen) to each other. Apart 
from the fact that this will establish a different proof of the Conley-Zehnder inequalities, 
this will also largely clarify the orbit structure of the auxiliary gradient-like system. 

Proof strategy. We will construct a particular homotopy of vector fields Vg, with s G [0, 1], 
that generate gradient-like flows. Their zeros will be independent of s, and hyperbolic. The 
vector fields interpolated by Vs are, for s = 1, fi = VgH, so that the zeros of Vg are precisely 
the critical points of H, and for s = 0, fo is a vector field that is tangent to every Cgen- For 
every s G [0, 1], we will construct a coboundary operator via the one- dimensional integral 
curves of Vg that connect its zeros. These coboundary operators are independent of s, and 
act on the free Z-module C of the Morse- Witten complex associated to (M, H). The desired 
estimate ( ^0[) then follows from a simple dimension argument. 



Construction of vq. We require vq to be gradient-like, and tangent to €gen- Furthermore, 
the zeros of Vo shall be hyperbolic, and shall coincide with the critical points of H. Con- 
sequently, the dimension of any unstable manifold of the flow generated by —Vq equals the 
Morse index of the critical point of H from which it emanates. 

To this end, we use the vector field constructed in the proof of proposition [11.6. 



We recall that it has been obtained by suitably deforming ny^gH such that all elements (ti 
of £ \ (tgen are removed. Furthermore, let us introduce compact e-neighborhoods Ue{^gen) of 
the generic connected components of C in the way demonstrated for ( p5|). 

It is now possible to extend the projector Q : T^^^^M —> T<Lgen that has been introduced 
in [11.2. 1| over the whole embedding space TUt{Cgen) ■ To this end, we pick an arbitrary smooth 
distribution W over the base manifold Ue{^gen), whose fibres over Cgen shall coincide with 
the corresponding fibres of TCgen- Admitting a slight abuse of notation, we let Q denote the 
^f-orthogonal projector TM W. Clearly, evaluating Q in any a G Cgen gives the projector 
Qa '■ TaM TaCgen discusscd in |11.2.1| . Because W is smooth, Q and its orthogonal 
complement Q are both smooth tensor fields. 

We define the vector field vq by requiring that it shall equal in M \ [/^(^gen), and 
that for X in U^{<tgen), it shall be given by 

vo{x) = {7rvVgH){x) + h{x){QVgH){x), 

where h : Ue{^gen) [0,1] is a smooth function obeying h\if^^^ = 1 and /i|ac/,(e;g^„) = 0. 
In particular, h shall be strictly monotonic along all non-constant trajectories of the flow 
generated by nyVgH, and the one form dh shall vanish on €gen- 

It can now be easily verified that Vq possesses all of the desired properties. That it 
generates a gradient-hke flow can be seen from the fact that outside of Ue{^gen), gC^gH, vq) = 
gCVgHjX^) is strictly positive, as has been shown in the proof of proposition ^l.6\ . Inside of 
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Ue{(tgen), one finds 

giVgH,vo) = \\nvV,H\\l + h\\QVgH\\l, 

due to the ^f-orthogonality both of iry and Q. The first term on the right hand side vanishes 
everywhere on C^en, but nowhere else in U^{<Lgen)- The second term reduces to ||(5Vg-ff||^ on 
^gen- Since evidently, QV gH\(f^^^ is the gradient field of the Morse function H\^^^^ : (tgen R- 
relative to the Riemannian metric on T€gen that is induced by g, its zeros are precisely the 
critical points of H on C^en, and there are no other zeros apart from those. This shows that 
the right hand side of the above expression is strictly positive except at the critical points of 
H, where everything vanishes. Because the scalar product of Vq with VgH is strictly positive 
except at the critical points of H, it is clear that —vq generates a gradient-like flow ipQ^t, so 
that H is strictly decreasing along all non-constant orbits. Furthermore, it is also clear from 
the given construction that vq is tangent to Cgen- 

Next, we address the proof that the zeros of —Vq are hyperbolic, and that the number 
of negative eigenvalues of the Jacobian matrix at any zero equals the corresponding Morse 
index of H. 

To this end, we pick a local chart at a critical point a of H. The Jacobian matrix of 
f at a in this chart is given by 

JsiCa{vo) = JaCaircvVgH) + JaCaiQVgH) 
= Ip^PaiDlHf + IqMDIH^ 

= {DlHf + {IpP^-lQM{DlHf- (31) 

To explain this result, let us first of all notice that there is no dependency on the function h 
because dh\(fg^^ is zero. Furthermore, {D'^HY is defined as the matrix [g''^Hjk\a] in the given 
chart, and Pa denotes the matrix of TTv{a). Linearizing the vector fields nyVgH and QVgH 
at a critical point of H, all terms involving first derivatives of H are zero. This explains the 
second line. The third line simply follows from Qa = l2n — Qa- 

Normal hyperbolicity follows from the invertibility of JaCa(fo), which is a consequence 
of lemma [11. 1| below. 

Definition of the homotopy of vector fields. Let 

Vs = sVgH + (1 - s)vo 

with s G [0, 1]. We claim that for arbitrary s, —Vg generates a gradient-like flow ipg^t, in a 
manner that H is strictly decreasing along all non-constant orbits. Furthermore, we claim 
that the zeros of Vs are hyperbolic fixed points of ips,t that do not depend on s. It then follows 
that the dimensions of the unstable manifolds are, for all s, given by the Morse indices of 
the critical points of H from which they emanate. 

To prove these claims, we consider the scalar product 

g{VgH,Vs) = s\\\/gH\\l + (1 - s)gi\/gH,vo). 

The first term on the right hand side is obviously everywhere positive except at the critical 
points of H, and the same fact has been proved above for the second term. Thus, H is 
strictly decreasing along all non-constant orbits of ips^t, which proves that it is gradient-like. 



24 



The Jacobian of Vg at a critical point of H is given by 

Jaca{vs) = s{DlHf + {1 - s)JaCa{vo) 

= {DlHf + {l-s){Ip^Pa~lQMiDlHf 

If we can prove that JaCa(fs) is invertible for all s G [0, 1], it follows that the number of 
negative eigenvalues is independent of s. To prove that this is indeed the case, we observe 
that because {D^HY is invertible, one merely has to show that (l2n + (l — s)(/p^ Pa — Jg^Qa)) 
is invertible. This in turn is satisfied if {Ip^Pa — IqaQa) has no eigenvalues in [1, oo) c R. 

Lemma II. 1 The spectrum of Ip^Pa — Ig^Qa has empty intersection with [1, oo). 

Proof. By arguing by contradiction. Since Ip^Pa — IqaQa is selfadjoint with respect to Qa, 
it is diagonalizable, and has a real spectrum. Let us assume that k G [1, oo) is an eigenvalue 
associated to the eigenvector w G T^M, so that 

{Ip^Pa-lQaQa)W = KW. (32) 

We multiply both sides of the equahty sign with Pa^QaQa from the left, and get the equation 
PalQaQalpaPaW = {1 — k) Palq^QaW . On the other hand, multiphcation from the left with Pa 
gives Palq^QaU! = (1 + K)PaW. Combining these two results, one obtains Palq^QalpaPaW = 
(1 — K^)PaW, which shows that PaW is an eigenvector of Palq^QalpaPa '■ Va ^ Va that belongs 
to the eigenvalue (1 — k^). We have proved earlier that normal hyperbolicity of €.gen implies 
that Palq^QalpaPa is invertible on V^. In addition, it is evident that (1 — k^) < 0. 
Let us assume that PaW ^ 0. Then, 

gaiw,PaIqaQaIPaPaW) = QaiQalPaPaW , Qalp^PaW) > 

follows from the gf-orthogonality of the projectors. On the other hand, we also have 

ga{w, Palq^QalPaPaW) = (1 - K^)ga{w,PaW) = (1 - K^) Q a{PaW , PaW) < 0, 



which is a contradiction. Hence, PaW = 0. In this case, ( ^2]) reduces to Iq^QaW = kw, 
and multiplication with QalpaPa from the left gives Qalp^PalqaQaW = nQalp^PaW = 0. 
Because QalpaPa^QaQa is invertible on Na^gen, this implies that QaW is zero. Therefore, w 
is not contained in the intersection of the images of Ip^Pa and Iq^Qa- However, being an 
eigenvector that solves ( ^H), it must be contained in this space, which is a contradiction. ■ 



The conclusion is that for all s G [0,1], the zeros of Vg are hyperbolic fixed points 
of ips^ti and that the dimensions of the corresponding unstable manifolds are given by the 
respective Morse indices of H . 

Since dtips^t = —'^s depends smoothly on s, hence ips.t is in s. Thus, s smoothly 
parametrizes a homotopy of stable and unstable manifolds emanating from the critical points 
of H, which belong to the gradient-like flow ips,t- 
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Because the fixed points of tljs,t are independent of s, and since the corresponding 
dimensions of the unstable manifolds coincide with the Morse indices of the critical points 
of if, we again consider the free Z-module 

C = (BpC^ 

that is generated by the critical points of H, and graded by their Morse indices. 

For every fixed s, we define a coboundary operator on C, using the fiow 4's,t- In fact, 
picking a pair of critical points of H with a relative Morse index 1, consider the unstable 
manifold of a, and the stable manifold Wj~^, of a' associated to 'ips,t, which are both 
smoothly parametrized by s. Since s parametrizes a homotopy of such manifolds, they 
naturally inherit an orientation from the one picked for s = 1 in the definition of the Morse- 
Witten complex for (M, H). 

Let T^E denote regular energy surface for an arbitrary energy value E between H{a) and 
H{a'). The intersection of W^^ with any regular energy level surfaces S^; of if is transverse, 
because H strictly decreases along all non-constant orbits generated by —Vs- 

H'Ee and W^^, fl S^; are oriented submanifolds of S^;, and smoothly parametrized 
by s. Hence, they define two homotopies of manifolds in J^e- Their intersection number, 
being a homotopy invariant, is independent of s, hence it equals the value obtained in case 
of s = 1. This implies that the coboundary operators obtained for arbitrary s are identical 
to the 5-operator of the Morse- Witten complex given for s = 1, since their respective matrix 
elements are equal. 

To finally prove ( we use the fact that all stable and unstable manifolds of ipo^t 
are, by construction, either confined to some or otherwise, that they connect critical 
points lying on different Cj's. This is because the stable and unstable manifolds of the flow 
generated by nyVgH only connect different connectivity components of €gen- 

Let us next consider pairs of critical points of H with a relative Morse index 1 that lie 
on the same component €i G ^gen, and the corresponding stable and unstable manifolds of 
ilJo,t which are contained in (ti. Since vo\^ simply is the projection of Vgii|e^ to T£j, these 
stable and unstable manifolds are precisely those which were used to deflne the Morse- Witten 
complex on (Cj,iij). 

Picking only the stable and unstable manifolds of ipQ^t contained in Cgen? we construct 
an operator 6 acting on C in the same manner in which the coboundary operator was deflned. 
It is again a coboundary operator, but now it is given by 

6i denotes the coboundary operator of the Morse- Witten complex associated to the pair 

Finally, we denote by Pj : C ^ Cj the projection of the free Z-module C generated 
by all critical points of H to the one generated by those critical points which lie in Cj. The 
above construction makes it evident that removing all integral lines of —Vq that connect 
critical points on different connectivity components of Cgen, one arrives at 6i = PiSPi, so 
that 

6 = ®iPi5P, 
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(note that 6 can be written as 5 = (BiSPi). Inclusion of the missing integral lines would yield 
S, as our homotopy argument has proved. This immediately makes the inequality 

dim(im5|cp) > dim(im5|cp) 
clear. We observe that this is precisely what is expressed in the inequality ( ^U|). ■ 
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Ill STABILITY CRITERIA FOR FIXED POINTS 



In this section, we will be interested in formulating stability criteria for equilibrium solutions 
of the constrained Hamiltonian system (M, u, H, V). In case of exponential (in) stability, the 
discussion is elementary, and the results are standard. 

However, in the marginally stable case, in which the linearized dynamics is oscillatory 
in nature, it is much harder to arrive at stability criteria if V is not integrable (the integrable 
case is again not interesting for us). 

We will give a heuristic line of arguments that supports a certain stability criterion 
for this situation. It will be obtained from an elementary application of averaging theory, 
and involves an incommensurability condition imposed on the frequencies defined by the 
linearized problem. Using a perturbation expansion that is adapted to the flag olV , we will 
argue why this condition, which could merely be an artefact of the averaging method, can 
presumably not be dropped. 

A rigorous proof of the conjectured stability criterion is far beyond the scope of this 
text, and is presumably at least as hard as proofs in KAM and Nekhoroshev theory. 



III.l Results of averaging 

For the discussion of stability, we again recall use auxiliary Kahler metric g on M, and denote 
the induced Riemannian distance function by distn (in contrast to the Carnot-Caratheodory 
distance function dc-c induced by g, which will be considered later). 

Definition III.l A point xq E € is stable if there exists 6{e) > for every e > 0, so that 
for all t, distFi{^t{x),xo) < e holds for all x with distR{x,xo) < S{e). 

Let xq e ^gen- and pick some small neighborhood U{a) C M together with an associ- 
ated Darboux chart, with its origin at Xq. The equations of motion are given by 

dtx = P{x)JH,{x) = X^{x), (33) 

where x = {x^, . . . ,x'^,Xn+i, ■ ■ ■ ,X2n), and JT" is the symplectic standard matrix. Further- 
more, Hx abbreviates d^H, and P is the matrix of the projector ny- a;-orthogonality of ny 
translates into P{x)J'X{x) — J'P'^{x)X{x) for all vector fields X. 

We will next bring the equations of motion in the vicinity of 0, corresponding to the 
point Xq, into standard form. To this end, let 

u := PqX y := PqX 

denote new coordinates, picked to be mutually orthogonal with respect to the Kahler metric 
5f|o in ^ TqM. 

Lemma III.l Locally, €gen is the graph of a function F : PqR'^'^ — > R^'^, u i— > F[u\. 
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Proof. The images of Pq and To€gen together span TqM = R^". To see this, let Qq again 
denote the orthoprojector TqM — > Tod^en- The claim is then implied by the fact that the 
matrix Ip^Po + IqqQq — l2n + {IpqPq ~ IqoQo) is invertible. The latter has already been 
proved above in the discussion of the induced Morse- Witten complexes on ^j. Thus, the rank 
of -Pq/qoQo = PoilpoPo + IqqQo) equals the rank of Pq, namely 2{n — k), hence Pq projects 
TqM surjectively onto To(tgen- Consequently, there is a linear map G : PqH'^'^ R^^ for 
which TqM is the graph {u,G{u)). (tgen thus admits a local parametrization {u,F[u]) with 
F{u) — G{u) + for sufficiently small \\u\\, which is since €gen is assumed to be 

smooth. ■ 



Proposition III.l The equations of motion are, in the coordinates {y,z), given by 

dty{t) = EQy{t) + Y[z{t),y{t)] (34) 
dtz{t) = Z[z{t),y{t)], (35) 

where \Y[z,y]\,\Z[z,y]\ ^ 0{\\y\\ \\z\\) + 0{\\y\\^), and Eq := PoM^[0, 0]7p„. 

Proof. We introduce the function y[y, u] := y — F[u], which defines the coordinate transfor- 
mation $ : {u,y) I— > {z,y). The Jacobi matrix of its inverse at {z,y) is 



D<^-'[z,y] = 



PQ + DF[u{z,y)]PQ 
-DF[u{z,y)]PQ Pq 



Prom the definition of F follows that PoDF[u\ — DF[u\. Thus, the equations of motion are 
now represented by 

dtz = (Po + DF[u{z,y)]PQ) X[z,y] (36) 

and 

dty = {Pq - DF[u{z,y)]Po) X[z,y] , (37) 

where = Xj^{^~^{z,y)). In this chart, {z,0) parametrizes Cgen, thus, by definition of 

<tgen: X[z, 0] = for all z. Taylor expansion of relative to (0, 0) gives 



= $-^[0,0] + £'$-^[0,0] ^ J +0(11x11^). 
Because of F(0) = 0, $~^(0, 0) = (0, 0). Taylor expansion of X^{u, y) relative to (0, 0) yields 

Xl[y,z] = M^[0,0] (^"t^ + R[u,y] , 
where R[u,y\ is a quadratic remainder term. Consequently, 

X[z,y] - M^[0,0] D$-i[0,0] ^ ^+R[y.z] 
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with a quadratic remainder term R[y, z\. It follows from Po-DX)^(0, 0) = that ( ^) holds 
for 

Z[z,y] = {Po + DF[u{z,y)]Po)R[z,y]. (38) 

Furthermore, 

dty = PoM]^[0,0] Jpoy 

-Vrn nl n „,M / ^ 



- Po M]^[0,0] DF[u{z,y)] ^ J + (39) 

where Y[z, y\ is a quadratic remainder terms. Finally, the kernel of DX)^[0, 0] is the tangent 
space ToCgen, which is also the image of DP[0], thus Z)X]^[0, 0]Z)P[0] = 0. Thus, 

dty = Po M]^[0,0] Ip, y + Y[z,yl (40) 

with a quadratic Taylor remainder term ^[2;, y]. Clearly, Z)X]^[0,0] = JIp^^Aq. 

^gen H U{xo) is a ccutcr manifold for the local system of ordinary differential equations 
( ^OD and ( |3^), which is parametrized by {z,y = 0). ■ 



Asymptotic (in) stability. A standard application of the center manifold theorem shows 
that if the spectrum of Hq does not intersect iH, there is a coordinate transformation {y, z) —>■ 
(y, z), so that ( ^) and ( ^8|) can be written as 

= + Y[yit),z{t)] 

dtz{t) = (41) 

PU| , where Y{0,z) = for all z. Thus, Xq G C^en is asymptotically unstable if there are 
eigenvalues with a positive real part, and asymptotically stable if all eigenvalues have a 
negative real part. 

In case of integrable V, asymptotic stability is impossible, because the eigenvalues al- 
ways come in pairs or quadruples with both positive and negative real parts. However, if 
V is nonintegrable, there is to the author's knowledge no obstruction to the existence of 
asymptotically stable equilibria, as the flow map is not symplectic. 



Marginal stability. The case of marginal stability is given when spec(So) C zR \ {0}. 

We will in this case conjecture a stability criterion, which is motivated by heuristic 
results obtained from averaging theory on the one hand, and from a perturbation expansion 
that is adapted to the flag of V, on the other hand. A rigorous stability analysis is beyond 
the scope of this text. 

To this end, let us assume that Sq is diagonalizable over C. To diagonalize it, we 
assume that the vector fields and Z[y,z] in ( ^) and ( ^) are analytic in {y,z), 

so that they possess a unique analytical continuation into a complex vicinity of xq G Cgen- 
The system of ordinary differential equations ( ^) and ( ^8]) is then interpreted in the sense 
that {y, z) is a vector in C'^^ x C^"~^^ (of small norm). The continuation of ^gen into C^" is 
defined by the common zeros of F [0, z\ and Z[0, z\ for z G C^*-"~'^\ 
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Let spec(So) = {iuji, . . . ,iuj2k}, with Ui G R. There exists a hnear transformation 

such that So is diagonal in the new coordinates, which, by abuse of notation, we denote 
again by {y, z). The equations of motion then reduce to 

dtyit) = di^g{iij)y{t) + Y[y{t),z{t)] 

dtz{t) = Z[yit),z{t)], (42) 

where u denotes the vector (cui, . . . , uj2k)- 

Let us next introduce polar coordinates (/, (p) and (J, 9) in terms of 

with r = 1, . . . , 2fc and s = 1, . . . , 2n-2A;. In particular, / e R^'', J G R^"-^'^, e [0, 27r]'^'' = 
T^^ (the 2A;-dimensional torus), and 9 e [0, 27r]^"-~^^ = T^"~^^. For brevity, vectors (e*'^''f'") 
and (e^^^w^) will be denoted by e^'^v and e*^iu, respectively, where v G R^'^ and if G R^"~^'^. 

In polar coordinates, the complexified equations of motion for y (the dot abbreviates 
dt) are given by 

e'^i + i^e^'^I = diagiiuj)e''^I + Y[e''^I, e'V], (43) 

so that 

/ = Re{e~'^Y[e'^I,e'''j]} 

= u; + Im{e-^'^diag(57)y[e^'^J,e^V]}. (44) 

In the same manner, 

9 = Im{e-^^diag(9j)Z[e*^J,e*''j]}. (45) 

(/, J) G R^" lies in a small vicinity of the origin. 

Next, we fix a small parameter e := ||/(0)||, and require that ||J(0)|| < O(e^). Then, 
we redefine the variables I tl and J — > t^J by rescaling. The new coordinates (/, J) have 
a norm of the order 0(1). 

Analyticity of Y[y, z\ and Z[y, z\ in (y, z) implies that the power series expansion with 
respect to e**^/ and e*^ J converges for all (/, J) G R^" sufficiently close to the origin. In this 
manner, ( 0) and ( ^5]) yield 

|m| + |p|>2 

js ^ J2 el'"l+2|^l-2G^p(/, J)e'<'"''^>e'<P'''> (47) 

|m| + |p|>2 

0, = a;,+ J2 el™l+2bl-i$^,.^p(j, j)e*(W)-</'-)e^(P'^) (48) 

|m| + |p|>2 

9, = J2 e''"'"^''^'~'e,;^p(J, J)e'<^'"''^>e'<P''\ (49) 

|m| + |p|>2 
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introducing the multiindices m G Z^'^ and p G Z^"^^^, with \m\ := ^ Im^l and \p\ := ^2 \Ps\- 
Of course, the right hand sides here are simply the Fourier expansions with respect to the 
27r-periodic angular variables and 6. Every Fourier coefficient, labelled by a pair of indices 
{■m,p), is a homogenous polynomial of degree \m\ in J, and of degree \p\ in J. 
In the limit e ^ 0, 

/ = j = ^ = uj 9 = . (50) 

Under the assumption that all components of uj are rationally independent, averaging can 
be applied with respect to the variable in order to obtain an approximation to the long 
time behaviour of the perturbed system. 

The new, "averaged" variables are obtained via 



If the frequencies Ur satisfy certain incommensurability conditions, it can be proved that the 
long time behaviour of the system is very well approximated by the averaged solutions. 

The only quantities in ( ^) ~ ( |4^) that do not vanish under averaging with respect 
to are the corresponding zero mode Fourier coefficients. 



We recall that in the initial equations of motion ( ^2D, the functions Y[y, z] and Z[y, z] 
are at least 0(||?/||). Thus, the leading terms of their power series in {y,z) are at least 
homogenous of degree 1 in y, and thus involve terms e^^'^''^'' with |m| > 1, but no terms with 
\m\ = 0. A quick inspection of the right hand sides of ( ^4]) and ( assures us of the fact 
that ( ^) and ( |4^) yield terms that survive the averaging process, but not ( and ( [491) . 



Therefore, averaging the perturbed equations of motion ( ^6]) ~ ( ^) with respect to gives 

I = e^F{I,J,e) J = ^ = (51) 

for some function F, where the bars account for averaged variables. 

Returning to the initial coordinate chart for real {y, z), and the notation xq for the fixed 
point in discussion, let us now assume that the quadratic form on Vx^ defined by Dl^H\v^^ 
is positive definite. We claim that the averaging result then suggests that xq is stable. To 
this end, we notice that Taylor expansion of the Hamiltonian (which is an integral of motion 
for the constrained system) relative to xq gives 

H{x) = H{xo) + u'dx.H{xo) + ^y'y' dx^dx,H{xo) + 0{e^), 

and recall that u = Px^x. Because of lemma [1II.1| the assumption < O(e^) implies that 
11^11 < O(e^). The quantity H{x) — H{xq) = O(e^) is an integral of motion; thus, if the 
quadratic form on Vx^ defined by D^^H is positive definite, \\y\\ has an order of magnitude 
0(e) as long as = O(e^) remains valid. The averaged equations of motion imply that 
II ^11 = II ^11 = ^(^^) is time- independent, so that due to -u = PqZ, the same applies to 
In conclusion, we conjecture the following stability criterion. 

Conjecture III.l Let xq G ^gen, spec(Ho) = {iuJi, . . . ,iuj2k}, with cjj G R \ {0}. Assume 
that (1) the frequencies uJr are rationally independent, and (2) that the quadratic form on 
VxQ defined by Dl^^H\y^^ is positive definite. Then, xq is stable. 
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III. 2 Dynamics along the flag of V 

Let us now approach the discussion of marginal stabihty from a different angle. In the sub- 
sequent paragraphs, we develop a geometrically invariant description of the local dynamics 
in the vicinity of € which is adapted to the flag of V. Our motivation is to comment on the 
origin of the incommensurability condition imposed on the frequencies involved in the above 
stability criteria, and to give arguments why it can not be dropped. 

Let us assume that x lies in a small open neighborhood f/ of xq G C, and that C^en : = 



C n f/ satisfies the genericity condition of theorem [ILl 



Proposition III. 2 Let €gen = € Cl U have the genericity property formulated in theorem 
{ll-lj . Then, there exists e > such that every point x & U with dii{x, ^gen) < ^ is given by 

x = expgY{xo) , |s| < e 

for some Y G T(V) with \\Y\\gj^j < 1? a;o G C^en (^Ws^ denotes the 1-parameter group of 
diffeomorphisms generated by Y , with expg F =id). 
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Proof. We pick a spanning family {Fj G V{y)}\ of V , with = 1. If for all Xq G 

Txo^gen coutalus uo subspace of then 

exPi(tiYi + ...+t2fc>'2fe)(^9en) H U 

is an open tubular neighborhood of C^en in U, for tj G (— e, e). Because the normal space 
Nxo^gen IS dual to the span of the 1-forms dFi at Xq, this condition is satisfied if and only if 
the matrix [dFj{Yi)] = [Yi{Yj{H))] is invertible everywhere on ^gen- According to proposition 
II. IL this condition is indeed fulfilled. ■ 



Due to proposition |1II.2| , there is an element Y G T(V) with < 1, so that 

X = ^'e(xo) 

for some < e ^ L Since Xq G Cgen, it is clear that under the flow generated by X]^, 
^±t{xo) = Xq, thus the solution of ( ^ belonging to the initial condition x is given by 

^^*(xo) := $t o ^,(xo) = [^t o o (xo) . 

In particular, \E'* is the 1-parameter group of diffeomorphisms with respect to the variable e 
that is generated by the pushforward vector field 

Ytix) := ^t*Y{x) = d^toYi^^tix)) , (52) 

where d^t denotes the tangent map associated to $t. This is a standard fact of differential 
geometry, cf. for instance p2|. From the group property Ys+t = ^s*Yt follows that 



dtYt = ds\,^,^s.Yt = [X^,Yt] (53) 

holds everywhere in U. 
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We recall from the beginning of the previous section that there exists a local spanning 



family {Yi e r(V)}J^^ for V that satisfies 



with J :- 



Ik 
-Ik 



Furthermore, defining 9i{-) :— u!{Yi,-), 



where J*-' are the components of — — J. Finally, 

= ny{X^) = - Y,{H) r Y, 

in the basis {Y^'^^-^. 

Proposition III. 3 Let f,Fi E C°^{U), where Fi := Yi{H), i = 1, ... ,2k, and assume that 
Fi{^!\{xQ)), f{^{{xQ)) are real analytic m e. For X,Y e T{TM), let 

OyX = [y,...,[y,x]] 

denote the r-fold iterated Lie derivative. Then, for sufficiently small e. 



r>Q 



(54) 



Proof. Clearly, 



= - F,{¥,{xo)) Y,{f){¥Xxo)) 

= - F,(**(a:o)) J'" i^U n) (/ o **)(xo) . 



Using the Lie series 



(55) 



(56) 



we arrive at the assertion. 



Proposition III. 4 Assume that 1^=0 £ r(y), and let be the given local spanning 

family ofV. Then, Cy^j £ ^{Yi}, where Vi is the i-th flag element ofV. 

Proof. Because oi^t* ■ ^{V) r{V), Yt is a section of V for all t if it is for t = 0. The 
claim immediately follows from the definition of the flag of F. ■ 
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Proposition |111.4| implies that there are functions a*(t, ■) G C°°{U), z = 1, . . . , 2/c, so 

that 

Yt{x) = a\t,x)Yi. (57) 
Their time evolution is governed by the following proposition. 

Proposition III. 5 Let Yt=o = a^Yi define the initial condition, and introduce the matrix 

n,:=[Yi{F,){x)r'] . 

Then, pointwise in x, 

a"'{t,x) = {expi-tn,))]" ai + F,{x)Ri"'it,x)al, (58) 

where ^ 

ETit^x) := ds{eM-{t~s)n.,))^uj{[Yu^,,Y,],Yn^ . 

Proof. The initial condition at t = is given by Yq = a^Yi, that is, by a*(0, x) = ag. Thus, 
by the definition of Yf in ( ^2p, one has Yt = $j ^.Yi, so that 

a\t,x) Yi = ai ^t*Yi . 

From uj(Yi,Yj) = Jij, Jik = -Jh and JimJ"^^ = Sl, 

Now, taking the t-derivative on both sides of the equality sign, one finds 
dta^{t,x) = -4u;([X]^,i>t,F,] , Ffe) J^'^ 

= - 4 ($t *r.)(F,)(x) P uj {Yi , n) .1^'^ 

- 4 F,{x) P J^™ u {\Yu *Y.^ , n) 
= - a\t,x)Yi{Fj){x) J^"^ 

- a^o F^{x) P J'^^ u {[Yi, , n) . 

Using the variation of constants formula pointwise in x, one arrives at the assertion. ■ 



III. 2.1 Leading order perturbation theory 

Let us next use the small parameter e for perturbation theory. We will only be interested 
in a heuristic argument that demonstrates why the incommensurability condition on the 



eigenfrequencies in conjecture III.l is presumably necessary. 



The simplified case that we will consider is defined by the following assumptions: 
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(1) Qx = ^, constant for all x in U. 

(2) spec(n)= {iui, . . . , iu!2k}, with Ur G R. 

(3) II II := sup^ I I ^ 7. 

Let us briefly comment on the generic properties of {a;^}. Writing Q = JA, we decom- 
pose the matrix A = [Yi{Yj{H)){xo)] into its symmetric and antisymmetric parts A^ and 
A_, respectively. A_ = [[Yj,Yi\{H){xo)]/2 vanishes if V is integrable, as one deduces from 
the fact that Xh\xo is a vector in Vj^ for all Xq G ^gen, and from the Frobenius condition. 
The linear system of ODE's a = JA^a in the space of a*(t)'s is Hamiltonian, hence the 
spectrum of JA+, if it is purely imaginary, consists of complex conjugate pairs of eigenvalues 
in m (here, we have introduced the notation a := (a^, . . . ,a^'^)). If JA_ admits a small 
relative norm bound with respect to JA+, pairs of complex conjugate eigenvalues ±iu!r will 
generically be deformed in a manner that they lose the property of being identical up to 
sign. Thus, generically, we may assume that all frequencies Ur are distinct from each other, 
and there are as many negative as positive ones. 

From ( ^), one infers 

Yt = al (exp(-tfi))} Y, + Y1 ^d^D 

because |-Fj(x)| = 0(|a;|) = 0(e), since Fj{xo) = 0. 
Thus, 

[Yt,X] = aiexp{-tQ)) [Yi,X] + ^ 0(e) [Y„X] + ^ 0(1) 

for all X G r(TM), and x G Ue{xo). Assuming that everything is sufficiently smooth, 
iterating the Lie bracket Cvt r times produces 

(flair (exp(-tn));- + 0(e)) ■•■]], 

\m=l / 

plus a series of terms with less than r nested Lie commutators, which contribute to higher 
order corrections (that is, 0(e^+^)) of terms indexed by r' < r in ( ^). 

Let us, for a discussion of leading order perturbation theory along each flag element of 
V, drop the terms of order 0(e), and assume that everything is sufficiently smooth so that 
our considerations hold if t < 0(e~^). 

For fixed r, let us consider the term 

F^i^iixo))^' (£;y.)(/ovl/*)(a;o), (59) 

which describes the contribution of ( |5^) along the r-th flag component of V, at least for 
sufficiently small t,. 

To this end, we have 

F.i^lixo)) = Yt{F,){xo) + 0(e2) , (60) 
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due to -Fj(xo) = 0. Therefore, 



F,(^*(xo)) X' = e exp(-tn)7 + 0{e') , (61) 

straightforward calculation shows. 

Collecting all results obtained so far, the terms with r nested commutators in ( l5^) are 



T ! 



exp(-t fi)| fi;. \ \[ (exp(-t fi));- j . . . , [F,,, 1^] ■ ■ ■]](/)(xo) 

+ 0(e^+2), 

as long as (iist/{(\E'*(xo), Xq) < 0(e). This implies that for / G C°°{U), 
fmxo)) ^ f{xo) + J2'-T- f «o exp(-s fi)| fij X 

r>0 

X (j[ air {expi-s n)yf^ . . . , [Y,^, Y^] • ■ ■](/)(xo) , (62) 

up to errors of higher order in e for every fixed r, as long as (izstj?(\E'*(xo), a;o) < 0{e). 

If / is picked as the i-th coordinate function x\ so that /(\E'*(xo)) = x\t), the quantity 
[Fj^, . . . , [Yi^., Yi] ■ ■ ■]{f){xo) is the z-th coordinate of the vector field defined by the brackets 
at Xq. Consequently, ( ^2|) is the component decomposition of x*(t) relative to the flag of V 
at Xq, to leading order in e. 

By assumption for the simplified model, Q has a purely imaginary spectrum. In this 
case, the operator exp(— si7) has a norm 1 for all s. Consequently, the integrand of ( |6^) 
is bounded for all s. It follows that if the integral should diverge, and become bigger than 
0(e), it will take a time 

to do so along the flag element Vr- The leading term of order 0(e), corresponding to r = 0, 
remains bounded for all t, but higher order corrections to it might diverge. 
Let us next write 

a(s) = exp(— s Q) % 

2k 

= 5Z e„ exp{-iujas) , (63) 

a=l 

where {e^} is a orthonormal eigenbasis of u with respect to the standard scalar product in 
C^^, and spec(Jl) = {iuja}- The amplitudes G C are determined by the initial condition 
a*(t = 0) = Oq, and will be assumed to be nonzero. By linear recombination of the vector 
fields Yi, one can set = 6i^a- Then, ( |6^) can be written as 

E h,^.,...,Ut)[Y^,,...AY^.,Yl]■■■]{f){x,), (64) 
r>0 l;ii,...,ir 
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where 

Ii;h,...,ir{t) ■■= I ds uiAi j expj-is jo;; + ^ „ 

"^0 \m=l J \ \ m=l 

r \ ~^ / r 

771=1 / \m=l 

X <^ exp -it i^i + / ^ 



m=l 

Evidently, the nested commutators vanish if all indices ii, . . . , ir, I have equal values. 

III. 2. 2 Conditions for instability 

Let us introduce the set 

which we endow with the norm 

:= sup . 

i{t)eo(^Ht) 

Furthermore, let 

\\A\\ := sup {\Ai\}, 

i=l,...,2k 

where Ai are C-valued amplitudes. 

We will now determine in which situations ||J^''^(t)|| diverges in the limit t oo. 
To this end, let 

21 := {uji, . . . ,uj2k} 

and denote the r-fold sumset by 

21^ := 21 + ■ • ■ + 21 , 

^ V ' 

r times 

defined by the set containing all sums of r elements of 21. 

For two sets of real numbers 21 and 03, we define their distance as 

c/(2t,<B) := ini {\ai-bj\ | G 2t , 6^ G <B } . 

Then, it follows from ( 0) that if rf(2l^, -21) > 0, 

< c/(2l^, -21)-^ \\Q\ 
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(the sum over frequencies Y7m=i^jm ( element of 21^, and can only equal —ui 

if d{^r, —21) = 0). However, if (i(2tr, — 2t) = 0, there is a tuple of indices {/; ii, . . . , v} such 
that 

r 

/z;n,...,vW = -tuJiAi II A,^ . (72) 

m=l 

This is simply obtained by letting the sum of frequencies in ( |65|) tend to zero. Thus, in this 



case, 



||jW(t)|| ~ t, (73) 

that is, a divergence linear in t as t ^ oo (of course, the validity of our leading order 
perturbation theory breaks down as t ^ -). Only if there are simultaneously positive and 
negative frequencies, d{^r, —21) = is possible, but due to the remark at the beginning of 



subsection |III.2.1| , this situation must generically assumed to be given. 

Let us consider some examples. The fact that ||J^°'*(t)|| is bounded for all t is trivial. 
On the next level, r = 1, we consider the component in the direction of the first flag element 
Vi = [V, V]. The condition for the emergence of a divergent solution is that (i(2l, —21) = 0. 
This is precisely given if there is a pair of frequencies ±uji of equal modulus, but opposite 
sign. For r = 2, assuming that (i(2l, —21) > 0, the condition d{^2, —21) = implies that 
there is a triple of frequencies such that tUj^ + = —Ui.^, ij G {1, . . . , 2k}. If this occurs, 
the solution will diverge in the direction of the second flag element, V2 = [V, [V,V]]. The 
discussion for r > 2 continues in the same manner. 

In conclusion, we have arrived at the following proposition. 

Proposition III.6 //rf(2l^, -21) = for some r, then p^''\t)\\ = 0{t) fort^ 00. 

The above results suggest that if the frequencies of the linearized problem do not satisfy 
the incommensurability condition "(i(2lr, —21) > for all r", the equilibrium xq is unstable. 
However, it has also become clear that the time required for an orbit to leave a Riemannian 
e-neighborhood Ue{xo) is extremely large. Assume that d{'Qir, — 2t) = for some r < r{V) 
(the degree of non-holonomy of V). Then, it takes a time on the order of O(^) to exit f/e(xo) 
along the flag element according to our leading order perturbative results (due to the 
factor in (ID). 

Although the perturbative results lose their validity already after t < O(-), it is always 
possible to continue solutions by picking new charts, centered around new basis points after 
times of order 0(-^), say. There is no doubt that since (i(2tr, —21) = holds for all basis 
points x'o G Ue{xo) (by assumption, Q is constant in Ue{xo)), the same divergence will be 
observed in each chart, and the prediction that f/e(xo) will indeed be exited after a time of 
order O(^) is realistic. During this time, the orbit will not drift away from Ue{xo) fl (tgen 
in the direction of V^g if Dl^H\v^^ is positive definite on \4g, as required in the conjectured 
stability criterion. This suggests that the incommensurability condition imposed on the 
frequencies of the linearized system can indeed not be dropped. 

We claim that although in the light of Riemannian geometry, divergences of this type 
are completely unspectacular up to times of order > O(-), they are very severe in the light 
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of the natural, internal geometry of the system, which is not the Riemann, but the Carnot- 
Caratheodory structure induced by g. 



III. 2. 3 Relation to sub-Riemannian geometry 

The structure exhibited in this discussion clearly shows that the constrained Hamiltonian sys- 
tem (M, uj,H,V) has many characteristics of systems usually encountered in sub-Riemannian 
geometry §,0 0. 

The natural metric structure that accounts for the particular structure of the present 
system is obtained from the Carnot-Caratheodory distance function distc-c induced by the 
Riemannian metric g. It assigns to a pair of points x,y & M the length of the shortest 
^-horizontal ^-geodesic. 

If V satisfies the Chow condition, distc-c{x,y) is finite for all x,y & M. This is the 
Rashevsky-Chow theorem |19|. In this case, the Carnot-Caratheodory e-ball 



B^~^{xo) := {x e M \ distc-c{x,Xo) < e } 

is open in M. 

If V does not satisfy the Chow condition, pairs of points that cannot be joined by 
horizontal (7j\/-geodesics are assigned a Carnot-Caratheodory distance oo. In this case, 
M is locally foliated into submanifolds Nx of dimension {2n — rankl^(y)), with A in some 
index set, which are integral manifolds of the (necessarily integrable) final element K-(y) of 
the flag of V {r{V) denotes the degree of nonholonomy of V^). On every A^;^, the distribution 
V\ := j^V satisfies the Chow condition, where j\ : Nx — M is the inclusion. Therefore, all 
points X, y G Nx have a finite distance with respect to the Carnot-Caratheodory metric in- 
duced by the Riemannian metric ixQM- Every leaf Nx is an invariant manifold of the flow 

The ball-box theorem. Assume that the spanning family {yi^}^2^^ of TM is suitably 
picked so that {Yi^} spans the flag element K- Let the gf-length of all F^^'s be 1. Then, we 
define the 'quenched' box 



Box£(a;) := < expj 




K e (-1,1) 



in Nxi where A is suitably picked so that x G A^a- Evidently, if V satisfies Chow's condition, 
A^A = According to the ball-box theorem |T^, there are constants C > c > 0, such 
that 

BoXce(a;) C fif"^(x) C Boxce(a;) . 

Therefore, Carnot-Caratheodory e-balls can be approximated by quenched boxes of Rieman- 
nian geometry. 



Instabilities in the light of C-C geometry. The above perturbative results imply that 
if there is some r < r(V), for which (i(2tr, — 21) = 0, the flow $t blows up the quenched 
boxes, and thus the Carnot-Caratheodory e-ball around Xq G CCgen, linearly in t, and along 
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the direction of Vr- In fact, B^"^'"\xq) is widened along Vr at a rate linear in t. After a 
time of the order O(^), the Carnot-Caratheodory e-ball containing the initial condition is 
increased to an extent that it can only be embedded into a Carnot-Caratheodory Oe(l)-ball. 
In this light, the instabilities in discussion are very severe. 

This type of instability has no counterpart in systems with integrable constraints. 
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IV NONHOLONOMIC MECHANICS 



We will in this section focus on nonholonomic mechanical systems, and their relationship 

to the constrained Hamiltonian systems considered previously. The discussion is restricted 
to Hamiltonian mechanical systems that are subjected to linear nonholonomic constraints 
(Pfaffian constraints). 

Hamiltonian mechanics. Let {Q, g, U) be a Hamiltonian mechanical system, where Q is a 
smooth Riemannian n-manifold with a C°° metric tensor g, and where U G C°°{Q) denotes 
the potential energy. No gyroscopic forces are taken into consideration. Let g* denote the 
induced Riemannian metric on the cotangent bundle T*Q. For a given X e r(TM), let 
Ox be the 1-form defined by OxiY) = g{X,Y) for all Y e T{TQ). It follows then that 
g{X, Y) = g*{9x, Oy) for all X,Y e r{TQ). 

The Kahler metric of the previous discussion, which was denoted by the same letter g, 
will not appear again in the sequel. Therefore, writing g for the Riemannian metric on Q 
should hopefully not give rise to any confusion. 

In a local triviahzation of T*Q, a point x e T*Q is represented by a tuple {q\Pj), 
where are coordinates on Q, and pk are fibre coordinates in T*Q, with i, j — 1, . . . ,n. The 
natural symplectic 2-form associated to T*Q, written in coordinates as 

Uq = dq' A dpi = —d9o , 

i 

is exact. 9o = Pidq^ is referred to as the symplectic 1-form. 
We will only consider Hamiltonians of the form 

H{q,p) = \ 9*Mp) + U{q) E C^(T*Q) , (74) 

that is, kinetic plus potential energy. In local bundle coordinates, the associated Hamiltonian 
vector field Xh is given by 

i 

The orbits of the associated Hamiltonian flow satisfy 

q" = dp^H{q,p) , pj = -dgjH{q,p) . (75) 

The superscript dot abbreviates dt, and will be used throughout the discussion. 

Let Ai denote the space of smooth curves 7 : / C R — T*Q, with / compact and 
connected, and let t denote a coordinate on R. The basis one form dt defines a measure on 
R. The action functional is defined by X : ^7 ^ R, 

^ J/* (^*^o - i/07) (76) 
- J^dt (j^Piitmt) - H{q{t),p{t))) , 
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with 7 = +Pj"9p.). Denoting the base point projection by 

7r:T*Q^Q , 



let 

c := (vr o 7) : / > Q 

denote the projection of 7 to Q. We assume that ||c(/)|| is sufficiently small so that solutions 
of ( ^) exist, which connect the end points c{dl). Among all curves 7 : J — T*Q with fixed 
projected endpoints c{dl), the ones that extremize X are physical orbits of the system. 

Linear non-holonomic constraints. Let us next discuss the inclusion of linear constraints 
on a given Hamiltonian mechanical system. This is achieved by adding a rank k distribution 
W over Q to the existing data, and by specifying a physical law, the Holder variational prin- 
ciple, that generates the correct physical flow on T*Q. The orbits of the resulting dynamical 
system have the property that their projections to Q are l^-horizontal. The physics of the 
Holder principle is, for instance, discussed in 0]. 

Orthoprojectors. There is a (^-symmetric tensor p : TQ TQ with 

Ker(p) = , p{X) = X VX G T{TQ) . 

It will be referred to as the gf-orthogonal projector associated to W. We note that its matrix 
is in Mat(n x n,R), of rank k. W-^ denotes the (^-orthogonal complement of W. The 
^i-orthogonal projector associated to W-^ will be denoted by p, so that 

p + p = id . 

Furthermore, the dual of W, denoted by W*, is defined as the image of W under the 
isomorphism g : TQ T*Q, and likewise for (W*)-^ := g o W*. The corresponding g*- 
orthogonal projectors are denoted by p^" and p^", respectively. We use this notation because 
the matrices of the latter are represented by the transposed matrices of p and p in every 
standard coordinate chart. We will use the same symbols for the projectors and their ma- 
trices. 



Dynamics of the constrained mechanical system. Next, we derive the equations of 
motion of the constrained mechanical system from the Holder variational principle. The 



use of the orthoprojectors p, p is inspired by |11|. For a closely related approach to the 



Lagrangian theory of constrained mechanical systems, cf. . 

Definition IV. 1 A projective W -horizontal curve in T*Q is an embedding 7 : / C R — 
T*Q whose image c = vr o 7 under the base point projection it : T*Q ^ Q is everywhere 
tangent to W . 

Let 7s : / ^ T*Q, with s G [0,1], be a smooth one parameter family of curves for 
which the end points Cs{dl) are independent of s (where := tt o 7^). 
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Definition IV. 2 A W -horizontal variation of a projective W -horizontal curve ^ is a smooth 
one parameter family 7^ : R ^ T*Q, with s G [0, 1], for which ^(vr o 7^) is tangent to W , 
and 7o = 7. 

Let 



6q'{t) := ds\,^Qq\s,t) , Spk{t) := ds\^^oPkis,t) . 

To any VT-horizontal variation 7^ of a l^-horizontal curve 70 witli fixed projections of the 
boundaries 

(7ro7,)(9J) = (7ro7o)(9/), (77) 
so that Sq^loj = 0, we associate the action functional 

^bs] = ^ ($^P.(s,t)g'(s,t) - H{q{s,t),p{s,t))'^ dt . 



Theorem IV. 1 (Holder's principle) A projective W -horizontal curve 70 : / T*Q corre- 
sponds to a physical motion of the constrained mechanical system if it extremizesl[^s\ among 
all W -horizontal variations 7^ that satisfy f\7li). 

Hence, if 

5X[7s] = + J^Yl ((^' - d,.H)Sq' - {cf + dp^H)5p,) = (78) 

for all VT-horizontal variations of 70 that satisfy = 0, then 70 is a physical orbit. 

Theorem IV. 2 In the given local bundle chart, the Euler- Lagrange equations of the Holder 
variational principle are the differential-algebraic relations 

q = p{q)dpH{q,p) (79) 
p\q)p = -p\q)d,H{q,p) (80) 
p{q)dpH{q,p) = 0. (81) 

Proof. The boundary term vanishes due to Sq^\gi = 0. 
For any fixed value of t, one can write Sq{t) as 

k 
a=l 

where Ya is a ^f-orthonormal family of vector fields over c(/) that spans Wc{i)- Furthermore, 
fa G C°°(c(J)) are test functions obeying the boundary condition fa{c{dl)) = 0. 

Since fa and 6p are arbitrary, the terms in ( |78D that are contracted with 6q, and the 
ones that are contracted with 6p vanish independently. In case of 6q, one finds 

^ dtfa {p + dgH)iY^ =0 
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for all test functions fa- Thus, {p + dqH)iY^ = for all a = l,...,k, or equivalently, 
p^{p + dqH) = 0, which proves ( 

Since 70 is W^-horizontal, p{q)q = 0, so the 5p-dependent term in (5X[7s] gives 

j^dt{q- pdpHY ip^p), + jdt {pdpHy {p^p)^ = . 

The components of 6p in the images of p^(g) and p^{q) can be varied independently. Thus, 
both terms on the second line must vanish separately, as a consequence of which one obtains 
( |7S[) and (|HT|). ■ 

Definition IV. 3 The smooth suhmanifold 

^Phys := \ mdpH{q,p) = 0} C T*Q 

defined by (\8l\) will be referred to as the physical leaf. 



It contains all physical orbits of the system, that is, all smooth path 7 : R ^ ^phys C 
T*Q that satisfy the differential-algebraic relations of theorem IV.2| . 



Theorem IV. 3 Let H be of the form (74)- Then, there exists a unique physical orbit 
7 : R+ ^ Mphys with 7(0) = x for every x e Mphys- 



Proof. Our strategy consists of proving that the differential relations ( ^) ~ ( ^T]) define a 
unique section X of TTlphys, where VJlphys is regarded as an embedded submanifold of T*Q. 
Solution curves 7 : ^phys of •yit) = X\^(^t-^ fulfill ( |7^) ~ ( |81|) for all initial conditions 

7(0) G Tlphys- The assertion thus follows from the existence and uniqueness theorem of 
ordinary differential equations. 

To this end, let us cover Tlphys with local bundle charts of T*Q with coordinates {q,p). 
In case of the Hamiltonian ( |7^), the defining relation ( |81|) for Tlphys reduces to 

p{q)M~\q)p = M-\q)p\q)p = 0, 

where M denotes the matrix of g (M is determined by the mass distribution of the system), 
and where one uses the (yf-orthogonality of p. This implies that ( ^ip is equivalent to the 
condition p^{q)p = 0. Hence, Tlphys is the common zero level set of the n component functions 

Consequently, every section 

X = v^{q,p)dqr + Ws{q,p)dp^ 

of T'^phys is annihilated by the 1-forms 

rf(pV)i = dqr{p^p), dq" + dpXp'^Pjidps 

for z = 1, . . . , n (of which only n — k are linearly independent), on '^Xftphys- This is, in vector 
notation, expressed by the condition 

= {v'-dqr)p^p + {wsdp^)p^p 

= {V''dqr)p^p + p^W , 
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which shows that the components f of X determine the projection p^w, so that knowing 
the components v and p'^w suffices to uniquely reconstruct X. Consequently, the right 
hand sides of ( ^) and ( determine a unique section X of TTlphys, so that every curve 
7 : Wlphys, with arbitrary 7(0) G Tlphys, that satisfies 

dait) = x(7(t)) 

automatically fulffils ( |79]) ~ ( ^). This proves the assertion. ■ 



Equilibria. The equilibria of the constrained Hamiltonian mechanical system on Wlphys are 
obtained from the condition g = and p = in ( |7S|) ~ ( ^T]) , whereupon one arrives at 

p = , p\q)d,U{q)=0. 

The critical set is thus given by 

<tQ:={qEQ\p\q)d,Uiq) = 0}. (82) 

An application of Sard's theorem fully analogous to the proof of theorem [11. 1| shows that 
generically, this is a piecewise smooth, n — fc-dimensional submanifold of Q, (recall that the 
rank of p{q) is k). 

Symmetries. Let G be a Lie group, and let 

^ : G ^ Difr(g) 

h tph , = id , 
be a group action on the configuration manifold so that the following holds: 

(1) Invariance of the Riemannian metrics: g oiph = 9 and g* oipf^ = g* for all h E G. 

(2) Invariance of the potential energy: U o ipf^ = U for all h E G. 

(3) Invariance of the distributions: iph^W = W and ip^W* = W* for all h EG. 

Then, we will say that the constrained Hamiltonian mechanical system {Q,g,U.,W) 
exhibits a G-symmetry. 

IV. 1 Construction of the auxiliary extension 

Let us finally merge the nonholonomic mechanical system into a constrained Hamiltonian 
systems of the type considered in the previous sections. 

To this end, we will introduce a set of generalized Dirac constraints over the symplectic 
manifold {T*Q,ujq) in the way presented in section |. They will define a suitable symplec- 
tic distribution V", in a manner that the constrained Hamiltonian system {T*Q.,ujq,H.,V), 
with H given by ( |7^, contains the constrained mechanical system as a dynamical subsys- 
tem. Thus, {T*Q,ujq,H,V) extends the mechanical system in the sense announced in the 
introduction. An early inspiration for this construction stems from 



We require the following properties to be satisfied by the auxiliary constrained Hamil- 
tonian system (T*Q, uq, H, V). 
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(i) VJlphys is an invariant manifold under the flow $t generated by ( 

(ii) All orbits with initial conditions x G Tlphys satisfy the Euler-Lagrange 
equations of the Holder principle. 

(iii) VJlphys is marginally stable under $t. 

(iv) The critical set € of $t is a vector bundle over €.q, hence equilibria of the con- 
strained mechanical system are obtained from equilibria of the extension by base 
point projection. 

(v) Symmetries of the constrained mechanical system extend to symmetries of $t. 

Let us briefly comment on (iii) ~ (v). (iii) is of importance for numerical simulations 
of the mechanical system, (iv) makes it easy to extract information about the behaviour of 
the mechanical system from solutions of the auxiliary system. Condition (v) allows to apply 
reduction theory to the auxiliary system, in order to reduce the constrained mechanical sys- 
tem by a group action, if present. The choice for V is by no means unique, and depending 
on the speciflc problem at hand, other conditions than (iii) ~ (v) might be more useful. 

Construction of V. Guided by the above requirements, we shall now construct V. 
To this end, we pick a smooth, 5f*-orthonormal family of 1-forms {C/}/=i with 



Cikiq) dq^ , 




(83) 



It is clear that // G C°°{T*Q). 



(I) To satisfy conditions (i) and (iii), we require that the level surfaces 



^,:={{q,p)\fi{q,p)=firJ = l 



,...,n-k}, 



(84) 



with fi := (/ii, . . . ,/i„_fc), are integral manifolds of K(y)- Here, r(y) denotes the 
degree of non-holonomy of V, and evidently, 9Jto = ^phys- 

Condition (iii) is satisfled because 



is an integral of motion for orbits of $4. Since L grows monotonically with 
increasing and attains its (degenerate) minimum of value zero on Tlphys, it 
is a Lyapunov function for DJiphys- Anything better than marginal stability is 
prohibited by energy conservation. 
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(II) To satisfy condition (ii), we demand that p{q)q = 0, or equivalently, that 

Ci{q) = , V/ = l,...,n-fc, (85) 

shall be satisfied along all orbits {q(t),p(t)) of ( |^), owing to ( |7Up. 

(Ill) If the constrained mechanical system exhibits a G-symmetry, characterized by a 
group action ip : G ^ Diff(Q) so that iph*W = W Wh E G, the local family of 
1-forms {d} can be picked in a manner that ip^Ci = (i is satisfied for all h E G in 
a vicinity of the unit element e. Consequently, the functions fi{q,p) = h*{ChP) 
and their level sets 971^ are invariant under the group action. 

The condition that ( are integral manifolds of Vr(v) ^ ^ implies that all sections 
of V are annihilated by the 1-forms dfj, for / = 1, . . . ,n — k. Furthermore, the condition 



( ^5]) requires V to be annihilated by the 1-forms 

6 := CiM dq' + dps (86) 

s 

that are obtained from lifting (j to T*{T*Q), with I = 1, . . . ,n — k. 
Proposition IV. 1 The distribution 



V 



f|kerd/,j fl l^n^^^^^) C T{T*Q) 



is symplectic. 



Proof. V is symplectic iff its symplectic complement V-^ is. With the given data, the latter 
condition is more convenient to check. V-^ is locally spanned by the vector fields (Yi, . . . , Y2k) 
obtained from 

uo{Yj,-) = U-) , ^iYi+k,-) = dfj{-), (87) 

where / = 1, . . . ,k, and Uq = —dpi A dq\ 

V-^ is symplectic if and only if D := [^^(l/, Yj)] has values in GLYi{2{n — k)). 

Remark. In the present notation, capital indices range from 1 to A; if they label 1-forms, 
and from 1 to 2k if they label vector fields. 
In local bundle coordinates, 

dfi = {d,.fi){q,p)dq' + Cniq)M'^{q)dp,, 

where Mij are the components of the metric tensor g on Q, as before. Let us introduce the 
functions ii^(g) := [Cjj('?)] andF(g,p) := [9qj/i<:(g,]9)], both with values in MatR,(n x (n — /c)), 
which we use to assemble 



K :-- 







: T*Q — > MatR(2(n - A;) x 2n) . 
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Any component vector v : T*Q —>■ R^"" that locally represents an element of T{V) satisfies 
Kv = 0. The symplectic structure ujq is locally represented by J, defined in ( ||). One can 
easily verify that the J-th row vector of the matrix KJ'~^ is the component vector of Yj. In 
conclusion, introducing the matrices 



G{q) : = 
S{q,p) : = 

one immediately arrives at 



E\q) M-\q) E{q) 

F\q,p) M-\q) E{q) - E\q) M~\q) F{q,p) , 



D 



KJK 



G 

-G S 



Since (i has been picked a 5f*-orthonormal family of 1-forms on Q, it is clear that G{q) = In-k- 
Thus, D is invertible. This proves that V-^ is symplectic. ■ 

Construction of the projectors. Next, we determine the matrix of the cuo-orthogonal 
projector ny, which is associated to V, in the present bundle chart. Again, it is more 
convenient to carry out the construction for its complement first. 

Proposition IV. 2 The matrix of the ujq- orthogonal projector vfy associated to (consid- 
ered as a tensor field that maps r(T(T*Q)) to itself, with kernel V) is given by 



P 
T pt 



in the local bundle chart {q,p). The matrix T = T{q,p) is defined in (WW- 



Proof. The proof of lemma can be used for this proof. The inverse of 



IS 



S 

In— fc 



" In— fc 





where we recall that G{q) = In-k- The I-th column vector of the matrix KJ' ^ is the 
component vector of Yi (we have required that {Yi, . . . l2{n-fc)} spans V-^). This implies that 

Lemma IV. 1 The matrix of the g-orthogonal projector p associated to W-^ in the present 
chart is given by 



p{q) = M-\q) E{q) E\q) 



(89) 



Proof. The construction presently carried out for yfy can also be applied to p. One only 
has to replace V-^ by to W-^, and ujq by the Riemannian metric g on Q. An easy calculation 
immediately produces the asserted formula. The matrix of p is subsequently obtained from 
P + P 



1. For more details, cf. 11 
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Introducing 

T(g,p) := E{q)F^iq,p)piq) - p\q) F{q,p) E\q) 
a straightforward calculation produces the asserted formula for Tfy. ■ 



(90) 



Corollary IV. 1 In the given bundle coordinates, the matrix of is 

^y= _T pt 



where T = T{q,p) is defined in (WW- 
Proof. This is obtained from ttv + ttv = '^2n- ■ 
The cuo-orthogonality of nv is represented by 



vry(x) J = J nl(x) 



in the given chart. 



Theorem IV. 4 Let H be as in ( yJf^. Then, the dynamical system 



p 

-pT -T 



(91) 



corresponding to the contrained Hamiltonian system {T*Q,ujq, H,V) is an extension of the 
constrained mechanical system {Q, g, U, W). 

Proof. By construction, DJlphys is an invariant manifold of the associated flow <l>t, hence ( |81|) 
is fulfilled for all orbits of ( |9l|) with initial conditions in DJiphys- 

The equation q = pdpH in ( ^1|) obviously is ( |7^) . 

Next, using the notation / := (/i, . . . , fn-k)\ 

1 = M~'p, 

and substituting ( ^Oj) for T{q,p), the equation for p in ( |9l|) becomes 

p = -p^ dgH - E F^ + p+F/ . 

Since are invariant manifolds of the flow $t generated by ( ^l]), dtfi{q{t),p{t)) 
vanishes along all orbits of ( ^l]), so that F'^q + E'^M~^p = 0. This implies that 



p = -pt dgH + EE^ M-'p + pt dg Q/V) 



(92) 



RecaUing that p = M'^EE^ from (|9|), and using the fact that / = on ^ytphys, one arrives 
at ( RO) by multiplication with pt from the left. ■ 
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Theorem IV. 5 The critical set of (\9^ corresponds to 

<t = {iq,p) \qE€Q;pe (W;)^} . 
It is a vector bundle over the base space (tq, with fibres given by those of (W*)-^. 

Proof. Let us start with equation ( ^). As has been stated above, the second term on the 
right hand side of the equahty sign equals p^(q')p, and moreover, from ( one concludes 
that 

/V = IIpVII', • 

The Hamiltonian ( ^) can be decomposed into 

1 _ 2 

H{q,p) = H{q,p^p) + - ||pV||^* , 

due to the 5f*-orthogonality of p"'' and p"'", so that ( ^2|) can be written as 

p = - P^dgH{q, p^p) + p^p . 
The equilibria of ( ^1]) are therefore determined by the conditions 
pt(g)p = , pt(g)9,/f(g,ptp) = 0. 
Because H depends quadratically on p^p, the second condition can be reduced to 

pt(g) d,U{q) = 

using the first condition. Comparing this with ( ^2|), the assertion follows. ■ 

In particular, this fact implies that every equilibrium {qo,Po) of the extension defines 
a unique equilibrium go on €q, which is simply obtained from base point projection. 

Extension of symmetries. Let us finally prove that {T*Q,ujq, H,V) extends the G- 
symmetry : G DiS{Q) of the constrained mechanical system [Q, g,U,W), if there 
is one. To this end, we recall that the 1-forms C/ satisfy ipKj for all h E G close to the unit. 
Via its pullback, ip induces the group action 

^ := ^* : Gx T*Q — > T*Q 

on T*Q. This group action is symplectic, that is, "^X^o = ujq for all h E G. For a proof, 
consider for instance 

The 1-forms ^7, defined in ( ^), satisfy "^X^i = and likewise, //ovi'/j = fj is satisfied 
for all h E G close to the unit. The definition of V in proposition [IV. 1| thus implies that 

^h*V = V 

is satisfied for all h E G. Due to the fact that uj and V are both G-invariant, ttv and vfy are 
also invariant under the G- action ^! . 



51 



The Hamiltonian H in ( [7^ ) is G-invariant under by assumption on the constrained 
Hamiltonian mechanical system. Thus, Xh fulfills ^h*XH = Xh for all h & G, which implies 
that Xj^ = 71v{Xh) is G-invariant. 



Stability of equilibria. To analyze the stability of a given equilibrium solution qo E €q, 
it is necessary to determine the spectrum of the linearization of X^ at Xq = (go, 0). 

A straightforward calculation much in the style of the discussion above shows that in 
the present bundle chart, 



DX^ixo) 





-pWlUp - R 



where 



[R 



■jk\ 



[d^.Uipr^ip)ld,s{p)l.] e MatR(nxn) 



(93) 



Furthermore, D^JJ is the matrix of second derivatives of U . The conjectural stability crite- 
rion formulated in the previous section can now straightforwardly be applied to DXYi{xq). 



IV. 2 The topology of the critical manifold 

Let us now discuss the global topology of the critical set of the constrained Hamiltonian 
mechanical system, defined by 



q e CQ;pt(g)p = 0;pt(g)p = ^/i/C/(g) 



Here, (i is an orthonormal spanning family of one forms for the rank n — k annihilator of 
the rank k distribution W, and 

<^Q = {qe Q\p\q)d,Uiq) = O} 

is the critical set of the physical system on Tlphys- We recall that generically, is a smooth 
n — /c-dimensional submanifold of Q. Evidently, C is the smooth rank n — k vector bundle 

over the base manifold Cg, whose fibres are given by those of the annihilator of W. 

The arguments and results demonstrated in section two can be straightforwardly ap- 
plied to the present problem. First of all, we claim that €q is normal hyperbolic with respect 
to the gradient-like fiow ipt generated by 

dtq{t) = -p{q{t))VgU{qit)), 

and that it contains all critical points of the Morse function U, but no other conditional 
extrema of apart from those (it is gradient-like because along all of its non-constant 
orbits, ^U{t) = —g{p'V gU, pV gU)\g(^t) < holds, as p is an orthoprojector with respect to 
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the Riemannian metric g on Q). The fact that this is true can be proved by substituting 
M ^ Q, H ^ U, P ^ a, 5'(Kahier) 9, and C ^ in section two, and by applying 
the arguments used there. Hence, letting /Xj denote the index of the connected component 
Cqi, defined as the dimension of its unstable manifold, the Conley-Zehnder inequalities ( |2^) 
imply that for a compact, closed Q, the topological formula 

J2 X^-^^'dirnHPiCQi) = J2 A^dimi/?'(Q) + (1 + A) Q(A) (94) 

i,p p 

holds, where Cgj are the connected components of €q. Here, denotes the p-th de Rham 
cohomology group with suitable coefficients, and Q{t) is a polynomial with non-negative 
integer coefficients. The argument using the Morse- Witten complexes associated to {Q, U) 
and to (Cq, U\i!Iq) to derive ( 0) can also be carried out in the manner explained in section 
2.4. 

Our next issue is to discuss the global topology of C Clearly, C is not a compact 
submanifold of T*Q, hence the Conley-Zehnder inequalities of the second section, which 
were derived for compact, closed, generic critical manifolds, do not apply. However, since 
both €. and T*Q are vector bundles of a particular type, one can nevertheless prove a result 
that is closely related to ( ^4|). Nevertheless, we claim that given the above stated properties 
of Cq, the generalized Conley-Zehnder inequalities 

^ X^+^^^duaHl^^Li) = XPdimHP{T*Q) + (1 + A)Q(A) (95) 

i,p p 

are valid. In this formula, the connected components €i of € are vector bundles whose base 
manifolds are the connected components Cgj of Cg, and the numbers /ij are the indices of 
with respect to tpf The polynomial Q{t) exhibits non- negative integer coefficients, and 
H* denotes the de Rham cohomology based on differential forms with compact supports. 

In fact, ( |95|) is a straightforward consequence of the circumstance that the base space 
of any vector bundle is a deformation retract the vector bundle. Hence, €.q, being the zero 
section of is a deformation retract of €, and likewise, Q is a deformation retract of T*Q. 
Since the de Rham cohomology groups are invariant under retraction, one infers the equality 

H^{(t,) ^ HP{€q,) m{T*Q) = HP{Q). 

Hence, formula ( ^5]) is equivalent to the assertion that 

J2 A^+'^Mimi/f (CqO = J2 X^dimHP{Q) + (1 + A)Q(A). (96) 

i,p p 

But this has just been proved by application of CZ theory to the flow ipt on Q. The weak 
Conley-Zehnder inequalities derived from this result are hence given by 
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where Bp is the p-th Betti number of Q. In particular, for the special value A = —1, one 
obtains 

i,p i 

where x denotes the Euler characteristic. 



IV. 3 Numerically determining the critical manifold 

Let us finally formulate another possible application of the results so far for technical pur- 
poses. 

Knowledge about the location of equilibria is very crucial for the design of a constrained 
multibody system. It is particularly desirable to know whether a chosen set of parameters 
and constraints has the effect that the critical sets are generic or not. 

Technical systems in engineering applications are often very large, so that their equi- 
libria can usually only be constructed with the help of numerical methods. The analysis 
that we have employed in chapters two and four, aimed at investigating global topological 
properties of generic critical manifolds, inspires the following strategy. 

We claim that if C/ is a Morse function, whose critical points are known, and if Q is com- 
pact and closed, it is possible to numerically construct all generic connectivity components 
of Cq. 

This is because generic components of are smooth, n — A;-dimensional submanifolds 
of Q containing all critical points of U, and no other critical points oiU\(TQ. This information 
can be exploited to find sufficiently many points on €q, so that suitable spline interpolation 
permits the approximate reconstruction of an entire connectivity component. To this end, 
one picks a vicinity of a critical point a of U, and uses a fixed point solver to determine 
neighboring zeros of Iplq)^/ gU{q)\'^ , which are elements of €q close to a. Iterating this 
procedure with the critical points found in this manner, pieces of (^q of arbitrary size can be 
determined. 

If all critical points of U are a priori known, one can proceed like this to construct 
all connectivity components of €q that contain critical points of U. In this case, one is 
guaranteed to have found all of the generic components of € if the numerically determined 
connectivity components are closed, compact, and contain all critical points of U. 

We remark that the determination of the critical points of a Morse function [/ : Q — > R 
is a difficult numerical task. Attempting to find critical points by simulating the gradient 
fiow generated by —VgU is presumably very time costly, because the critical points define 
a thin set in M. Their existence, however, is ensured by the topology of Q if the latter is 
nontrivial, as a result from the Morse inequalities. 

Another remark is that all critical points a at which JaCa(pV(,t/) has a reduced rank, 
are elements of the non-generic part of €q. If there are such exceptional critical points in a 
technically relevant region of Q, they can be removed by a small local modification of the 
system parameters and the constraints. 
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